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ABSTRACT 


To  test  the  applicability  of  various  one-dimensional  sediment - 
transport  numerical  models  of  movable-bed  channels  to  Pool  20  of  the  Missis¬ 
sippi  River,  detailed  field  data  were  first  collected,  and  several  numerical 
models  were  subsequently  evaluated  using  these  field  data. 

Chapter  I  describes  the  detailed  results  from  a  sediment  sampling 
field  study  conducted  in  1978  near  Buzzard  Island  (EM  347-55)  in  Pool  20  of 
the  Mississippi  River.  Variations  in  the  longitudinal  and  transverse  dis¬ 
tributions  are  presented  for  various  flow  quantities  and  sediment  character¬ 
istics,  and  are  analyzed  in  connection  with  the  prevalent  shoaling  problem 
in  the  Buzzard  Island  study  reach. 

Chapter  II  evaluates  the  performance  of  four  mathematical  simula¬ 
tion  models  which  used  the  data  from  a  complimentary  field  study  conducted 
in  1976  to  establish  the  initial  conditions  when  predicting  the  1978  condi¬ 
tions.  The  four  models  evaluated  were:  (1)  the  HEC-6  model  provided  by  the 
Hydrologic  Engineering  Center,  U.S.  Army  Corps  of  Engineers;  (2)  the  SUSR 
and  (3)  L1IVSR  models  employed  at  the  Colorado  State  University;  and  (4)  the 
CIIAR2  model  used  by  the  Sogreah  consulting  firm  in  Grenoble,  France. 
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FIELD  STUDY  AND  TESTS  OF  SEVERAL  ONE-DIMENSIONAL 
SEDIMENT-TRANSPOTT  CXMPUTER  MODELS  FOR 
POOL  20,  MISSISSIPPI  RIVER 


1 .  1978  FIELD  STUDY  AND  DATA  PRESENTATION 


A .  Introduction .  The  field  study  was  conducted  to  obtain 
detailed  information  concerning  the  flow  and  sediment -transport  charac¬ 
teristics  along  the  Mississippi  River  reach  near  Buzzard  Island  (river 
mile  (RM)  349-51).  The  location  of  the  study  reach  is  shown  in  figure 
I-A.l.  The  field  data  collected  during  the  study  were  obtained  to  help 
explain  and  lead  to  a  better  understanding  of  the  river  processes  that 
create  a  shoaling  problem  in  the  vicinity  of  the  downstream  tip  of  Buz¬ 
zard  Island.  There  are  several  related  adverse  consequences  of  this 
shoaling  problem.  Foremost  is  the  adverse  effect  that  the  shoaling  has 
on  the  navigation  channel.  Significant  rises  in  the  bed  elevation  due 
to  sediment  deposition  can  obstruct  barge  traffic  within  the  navigation 
channel,  which  must  then  be  dredged  to  maintain  the  required  9-ft  depth 
along  the  passageway.  The  cost  of  dredging  has  risen  sharply  in  recent 
years,  making  dredging  an  expensive  means  of  alleviating  the  shoaling 
problem.  Another  problem  that  must  then  be  faced  pertains  to  the  disposal 
of  the  dredged  material.  This  dredged  material  must  be  disposed  of  in 
accordance  with  the  standards  established  by  the  various  environmental 
and  conservation  agencies,  so  as  not  to  be  detrimental  to  the  surrounding 
environment  or  wildlife. 

There  are  three  general  features  of  the  Mississippi  River  (MR) 
in  the  study  reach  that  lead  to  the  shoaling  problem.  First,  the  flow  bi¬ 
furcates  at  two  locations  near  the  downstream  tip  of  Buzzard  Island.  The 
bifurcation  of  flow  reduces  the  main  channel  flow  and  diminishes  the 
sediment-transport  capacity  of  the  river.  Second,  the  river  widens  in 
the  shoaling  area,  which  decreases  the  mean  flow  velocity  and  diminishes 
the  sedimnit-tninsport  capacity.  Finally,  the  thalweg  crosses  the  river 
in  t  his  reach.  This  so-eal  hxi  "cross-over"  reach  lacks  the  strong  secon¬ 
dary  currents  that  are  produced  in  river  bends  which  significantly 
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increase  the  sediment-transport  capacity  of  the  flow.  The  sediment 
being  deposited  in  the  study  reach  is  thought  to  originate  primarily 
from  the  Des  Moines  River  (DMR),  whose  confluence  with  the  MR  is  roughly 
10  mi  upstream  from  Buzzard  Island.  The  smaller  energy  slope  and  veloc¬ 
ity  of  the  MR,  in  comparison  to  that  of  the  DMR,  are  unable  to  carry  the 
sediment  load  that  the  DMR  transports  to  the  MR,  and  sediment  deposition 
occurs  downstream. 

The  three  main  objectives  of  this  study  were  as  follows.  The 
primary  objective  was  to  collect  data  on  the  transverse  and  streamwise 
variations  of  the  flow  depth,  flow  velocity,  suspended-sediment  and  bed¬ 
load  discharges,  and  bed-material  properties.  This  information  could  be 
used  to  analyze  the  previously  metioned  shoaling  problem.  Hie  second 
objective  was  to  employ  these  data  to  formulate  power-law  relation¬ 
ships  between  the  water  and  sediment  discharges.  These  relationships  were 
then  applied  to  estimate  the  change  expected  in  the  sediment  discharge  as 
a  result  of  any  change  in  the  flow  discharge,  for  example,  the  effect  of 
the  implementation  of  any  corrective  measure  designed  to  relieve  the  shoaling 
problem.  The  final  objective  was  to  obtain  complete  and  reliable  field 
data  which  could  be  utilized  in  the  formation  and  testing  of  mathematical 
sediment-transport  simulation  models. 

The  final  objective,  regarding  the  evaluation  of  numerical  simu¬ 
lation  models,  is  the  topic  presented  in  Chapter  II.  It  should  be  mentioned 
that  the  1978  field  study  near  Buzzard  Island  was  done  to  complement  an 
earlier  field  study  conducted  in  1976.  The  earlier  study  (Nakato  and 
Kennedy,  1977)  investigated  the  sediment  transport  characteristics  of  the 
MR  near  Fox  Island  (RM  355-6)  and  Buzzard  Island  (RM  449-50).  As  will  be 
discussed  in  Chapter  II,  the  data  collected  during  the  1976  field  study 
were  used  to  establish  the  initial  conditions  for  four  numerical  models. 

The  models  were  then  run  for  the  28-month  time  period  between  the  1976  and 
1978  field  studies,  and  the  model  predictions  of  the  river  characteristics 
were  compared  to  results  found  during  the  latter  study. 

The  text  that  follows  presents  a  condensed  version  of  the  origi¬ 
nal  report  that  was  prepared  for  the  1978  field  study  (Vadnal ,  1978). 

Readers  who  desire  further  clarification  or  more  detail  concerning  the 
various  subjects  discussed  lx 're  in  should  refer  to  the  original  ro|*>rl.. 


B.  Equipment  and  Field-Laboratory  Procedure.  Data  collection 
on  the  MR  was  conducted  from  the  Iowa  Institute  of  Hydraulic  Research 
(IIHR)  18-foot  Jon-Boat.  A  Price  AA  bucket -type  current  meter  was  used 
to  obtain  flow  velocities,  while  the  bed-material,  depth-integrated 
suspended-sediment,  and  bed-load  samples  were  collected  using  a  US  BM-54, 

US  P-61,  and  Helley-Smith  sampler,  respectively.  Sediment  samples  were 
taken  at  eight  verticals  that  were  fairly  evenly  spaced  across  the  main- 
channel  cross  sections  of  the  MR,  while  typically  3  or  4  verticals  were 
sufficient  for  sampling  the  side-channel  sections  and  the  DMR.  Consequently, 
these  verticals  established  an  identical  number  of  subsections  at  each 
cross  section  for  which  the  flow  and  sediment  discharge  characteristics 
would  be  determined.  In  addition  to  the  sediment  samples,  flow  velocities 
were  measured  at  each  vertical  at  two-tenths  and  eight -tenths  of  the  flow 
depth,  which  were  then  averaged  to  obtain  a  mean  subsection  flow  velocity. 
Detailed  measurements  of  flow  velocities  and  suspended-sediment  concentra¬ 
tions  were  made  at  one  vertical  for  each  main-channel  section  in  the  MR 
and  for  the  DMR. 

Analyses  of  all  sediment  samples  were  performed  at  the  IIHR 
according  to  reoomnended  laboratory  procedures  (Guy,  1969).  The  analyses 
of  the  suspended-sediment  samples  gave  the  suspended-sediment  concentrations, 
in  parts  per  million  (mg/1).  Analyses  of  the  bed-load  samples  gave  the 
bed-load  discharges  and  the  bed- load  particle-size  graduations,  while  the 
bed-material  analyses  gave  the  bed-material  particle-size  distributions. 

The  results  of  the  particle-size  analyses  of  all  the  bed-material  samples 
collected  during  the  two  trips  of  the  1978  field  study  are  tabulated  in 
Appendix  A. 


C.  Presentation  and  Discussion  of  Field  Data. 

1.  Data  analysis.  Locations  of  twelve  sampling  sections  were 
chosen  in  an  arrangement  that  would  enable  the  determination  of  the  flow 
and  sediment  characteristics  of  the  MR  in  the  Buzzard  Island  study  reach. 
r[\w>  dnta-eol looting  trips  which  gathered  a  total  of  961  sediment  s;unplos 
were  nude  during  (he  I97S  field  study:  Mm*  first  trip  during  2C>  June  - 
112  July  ;md  Llie  s<vond  trip  during  14-25  Aujru^t..  'the  total  disc!  large  at 
each  cross  section,  Q,  was  determined  by  sunvning  the  products  of  each 
average  subsection  flow  velocity  and  the  subsection  area,  determined 
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from  the  cross-section  plot.  Similar  procedures  were  used  to  compute 
the  suspended-sediment  discharge,  Qg,  and  bed-load  discharge,  Q^.  The 
resulting  size  distributions  of  the  bed-material  samples  were  plotted  on 
semi-logarithmic  graphs  to  evaluate  the  median  diameter,  D^g,  geometric 
mean  diameter,  (=  /Dg^D-^g),  and  geometric  standard  deviation, 

(=  /Dg^/D^g).  The  value  of  D^g  represents  the  particle  size  for  which 
50  percent  of  the  sample  by  weight  is  finer;  similar  definitions  apply 
to  and  Dl6. 

A  sunmary  of  the  principal  quantities  obtained  from  the  field 
data  is  presented  in  table  I-C-l.l.  River  cross  sections  are  labelled 
with  three  digits.  The  first  refers  to  the  physical  location  of  the 
cross  section.  The  second  digit  locates  the  cross  section  with  regard 
to  passage  around  an  island.  The  channel  closer  to  the  right  river  bank 
is  (1)  and  the  one  closer  to  the  left  is  (2).  The  third  digit  describes 
whether  data  were  measured  on  the  first  or  second  trip.  The  first  and 
second  sampling  trips  are  identified  by  1  and  2,  respectively,  as  the 
third  digit  in  sequence  and  are  listed  on  separate  pages.  Other  para¬ 
meters  listed  include  water  temperature,  T,  energy  slope,  S  (determined 
from  the  upstream  and  downstream  pool  elevations  recorded  at  Lock  and 
Dams  (L&D)  19  and  20,  respectively),  top  river  width,  W,  mean  flow  depth, 
d(=A/W),  and  mean  flow  velocity  U(=Q/A),  where  A  is  total  cross-section 
flow  area. 


2.  Cross-section  profiles.  Cross-sectional  profiles  of  the 
main-  and  side-channel  sections  compiled  from  the  first  and  second  trips  are 
shown  in  figure  I-C-2.1.  Five  cross  sections,  which  will  be  used  later 
to  evaluate  mathematical  model  performances,  were  selected  out  of  the  12 
cross  sections  for  their  proximity  to  the  1976  field-study  cross  sections. 
These  cross  sections  were  used  to  fix  the  initial  conditions  in  the  models. 
Table  I-C-2.1  lists  the  locations  of  the  cross  sections  from  the  1976  and 
1978  field  studies  that  will  serve  as  the  bases  for  comparisons  in  Chapter 
II. 


Table  I-C-l.l  Summary  of  Principal  Quantities  Obtained  from  the  Field  Data (1st  Trip) 


Table  I-C-l.l  (cont*d)  (2nd  Tri 
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Table  I-C-2.1 


Comparable  Cross  Sections 

from  the  1976  and  1978 

Field  Studies 

1978 

BM 

1976 

RM 

1-2 

354.89 

3-2 

354.93 

6-1 

349.75 

5-1 

349.82 

7-1 

349.50 

6-1 

349.45 

8-2 

349.24 

7-2 

349.29 

9-1 

348.98 

8-1 

348.96 

3.  Lateral  distributions.  Figures  I-C-3. 1  through  I-C-3. 26 
present  the  lateral  (cross-channel)  variations  in  depth,  d;  mean  flow 
velocity,  u;  unit  water  discharge,  q;  mean  suspended-sediment  concentra¬ 
tion,  C;  unit  suspended- sediment  discharge,  qg;  unit  bed-load  discharge, 
q^;  and  median  bed-material  diameter,  D  q,  obtained  from  the  field  data 
collected  during  both  trips  at  the  main- channel  sections.  Note  again 
that  when  referring  to  any  specific  section  (e.g. ,  1-2-1),  the  third 
number  in  the  sequence  denotes  the  trip  number.  All  cross  sections 
are  plotted  so  that  the  right  bank  is  on  the  left  side  of  the  figure,  as 
if  being  viewed  from  downstream. 

The  lateral  distributions  determined  for  the  first  trip  are 
shown  in  figures  I-C-3. 1  through  I-C-3. 13.  Figure  I-C-3. 1  of  section 
1-2-1  shows  a  high  mean  concentration  of  suspended  sediment  near  the  right 
bank.  This  high  concentration  was  caused  by  the  confluence  of  the  DMR 
and  the  MR  roughly  6  mi  upstream  from  this  section.  The  large  MR  flow 
velocities  swept  the  inflow  from  the  DMR  downstream  without  causing  much 
cross-stream  mixing.  The  DMR  inflow  carried  a  high  concentration  of  sus¬ 
pended  sediment  (approximately  1,800  ppm  when  measured  earlier  on  the  same 
day).  The  river  reach  between  the  confluence  and  section  1-2-1  is  fairly 
straight,  so  that  lateral  mixing  due  to  curvature  effects  is  minimal.  The 
figure  shows  that  the  concentration  of  about  60C  ppm  near  the  right  bank 
drops  across  the  channel  to  a  nearly  constant  concentration  of  about  200 
ppm.  Note  that  the  concentration  of  roughly  600  ppm  measured  near  the 
right  bank  was  one-third  of  the  concentration  measured  in  the  DMR.  The 
concentration  became  almost  uniform  across  the  channel  at  the  downstream 
section  8-2-1 ,  as  seen  in  figure  [-C-3.3. 
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The  lateral  distributions  at  section  2-1-1  are  shown  in  figure 
I-C-3.2.  This  sectioi  was  located  downstream  from  a  bend  in  the  river. 

The  centrifugal  effect  of  the  bend  and  secondary  currents  have  carved  a 
steep  right  bank.  Note  that  the  mean  flow  velocity  near  this  bank  for  the 
first  trip  was  almost  5  ft/s.  Tno  slight  increases  in  depth  and  mean  veloc¬ 
ity  near  the  left  bank  were  due  to  the  effect  of  the  flow  behind  the 
small  island  downstream  from  Taylor  Island.  The  large  drop  of  mean 
concentration  seen  in  the  shallow  center  of  the  section  was  due  to  the 
bifurcation  of  flow  downstream  at  section  3,  where  roughly  12  percent 
of  the  total  suspended-sediment  load  occurred  in  the  side  channel  section 
3-1-1. 

The  lateral  distributions  at  section  3-2-1  are  shown  in  figure 
I-C-3.3.  The  right  bank  of  this  section  is  located  on  the  upstream  tip  of 
Buzzard  Island.  Large  flow  velocities,  over  4  ft/s  near  the  right  bank, 
produced  a  strong  scouring  action,  making  the  tip  unstable  and  washing 
large  concentrations  of  sediment  downstream.  The  cross-sectional  profile 
of  section  3-2  in  figure  I-C-2.1  shows  that  a  lot  of  scouring  had  occurred 
near  the  right  bank  of  the  island  between  the  first  and  second  trips. 

During  both  trips,  the  tip  of  the  island  was  wedge-shaped  and  had  very 
steep  banks. 

The  lateral  distributions  at  section  4-1-1  in  figure  I-C-3.4 
show  extremely  high  suspended-sediment  concentrations  near  the  right  bank, 
some  of  which  were  caused  by  the  scouring  action  near  the  tip  of  Buzzard 
Island.  However,  the  extent  of  the  lateral  dispersion  of  suspended 
sediment  (seen  to  occur  over  1000  ft  in  a  longitudinal  distance  of  1900 
ft  between  sections  3-2  and  4-1)  seems  unrealistic,  so  another  sediment 
.source,  in  addition  to  the  scouring  of  the  island,  is  presumed.  Conse¬ 
quently,  the  major  cause  of  the  high  sediment  concentrations  was  most  likely 
the  suspended-sediment  inflow  from  the  DMR,  approximately  10  mi  upstream. 

A  mean  suspended-sediment  concentration  of  3,300  ppm  for  the  DMR  was 
measured  approximately  one  hour  after  the  sampling  at  section  4-1-1  was 
completed.  Note  that  the  mean  suspended-sediment  concentrations  of  roughly 
1,100  ppm  measured  near  fix'  right  bank  oT  section  4-1-1  are  one-third  of 
the  sediment  concentration  measured  in  the'  DMR.  Note  also  the  similar 
trend  described  for  section  1-2-1.  Although  the  unit  flow  discharge  was 
high  near  the  left  bank,  the  unit  suspended-sediment  discharge  was  much 
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larger  near  the  right  bank  where  the  high  mean  concentrations  were 
measured.  These  higher  concentrations  near  the  right  bank  continued  for 
approximately  two  miles  downstream,  and  finally  became  uniformly  distri¬ 
buted  across  the  transect  at  section  8-2-1.  Note  that  the  effect  of  the 
DMR  sediment  outflow  was  observed  to  extend  roughly  12  mi  downstream 
from  the  confluence. 

Another  pattern  that  developed  in  the  first  five  sections  was 
that  the  lateral  distributions  of  unit  bed  load,  q^,  were  all  single- 
peaked,  and  the  maximum  values  followed  the  navigation  channel  of  the  MR. 

The  distributions  became  multi-peaked  farther  downstream,  where  the  flow 
bifurcated  at  transects  7  and  9. 

The  lateral  distributions  at  sections  5  through  8  are  shown  in 
figures  I-C-3.5  through  I-C-3.8.  At  these  locations,  the  deepest 
parts  of  the  cross  sections  consistently  occurred  near  the  left  bank, 
away  from  the  navigation  channel  which  crosses  the  river  to  the  right 
bank.  The  major  bed-load  movement  at  section  6-1-1,  as  shown  in  figure 
I-C-3.6,  occurred  in  the  right  half  of  the  section,  especially  in  the 
region  just  beyond  the  submerged  wing  dam  that  extended  from  the  left  bank. 
The  bed-load  movement  appeared  restricted  in  the  left  half  of  the  section 
while  the  flow  velocities  and  unit  flow  discharges  decreased  quickly 
towards  the  left  bank.  Note  the  increase  in  the  magnitude  of  the  bed-load 
discharge  from  section  6-1-1  to  7-1-1  shown  in  figures  I-C-3.6  and  I-C-3.7, 
respectively,  with  a  similar  decrease  in  magnitude  from  section  7-1-1  to 
8-2-1  in  figures  I-C-3.7  and  I-C-3.8.  Sections  6-1-1  and  7-1-1  were 
sairpled  on  the  same  day,  while  section  8-2-1  was  sampled  on  the  following 
day.  The  reach  between  sections  6-1  and  7-1  will  be  discussed  later  con¬ 
cerning  scouring,  and  the  reach  between  sections  7-1  and  8-2  will  be 
discussed  with  shoaling. 

The  bifurcation  of  flow  occurring  at  transect  9  is  apparent 
from  figures  I-C-3.7  and  I-C-3.8  for  sections  7-1-1  and  8-2-1,  respectively. 
The  largest  values  of  depth,  mean  flow  velocity,  unit  water  discharge, 
and  unit  suspended- sediment  discharge  all  occurred  close  to  the  left 
bank.  Note  that  the  navigation  channel  is  aiong  the  right  bank.  Figure 
T-C-3.9  shows  the  lateral  distribution  at  section  9-1-1,  when*  the*  deepest 
part  ot  the  transect  is  near  the  right  I  Kink  when*  the  navigation  channel 
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is  located.  Table  I-C-3.1  summarizes  the  percentages  of  water  discharges 
and  sediment  discharges  through  side  channels  7-2-1  and  9-2-1. 


Table  I-C-3.1 

Water  and  Sediment  Discharges  at  Sections  7  and  9 


Section 

Date 

Q 

(cfs) 

Qs 

(tons/day) 

% 

(tons/day) 

7-1-1 

070578 

101,830 

89,467 

1,958 

7-2-1 

070678 

16,155 

13,280 

102 

Qs2/Qsi=15^ 

Q^/Qg!^ 

9-1-1 

070778 

96,800 

62,925 

630 

9-2-1 

070778 

30,271 

22,956 

77 

q2/q1=31% 

Qs2/Qsl=36% 

W125 

Figures  I -€-3.9  and  I-C-3.10,  which  depict  sections  9-1-1  and 
10-1-1,  respectively,  show  that  the  maximum  unit  bed-load  discharges  occurred 
near  the  left  bank.  The  flow  was  very  turbulent  near  the  upstream  tip 
of  Huff  Island  during  the  first  trip.  The  splitting  of  the  flow  around 
Huff  Island  and  the  presence  of  a  submerged  wing  dam  extending  from  the 
left  bank  upstream  from  section  9-1  created  turbulence  which  caused  the 
large  bed-load  movement .  The  figure  of  section  10-1-1  shows  that  the  mean 
velocity  was  fairly  uniform  across  the  transect ,  although  the  depth  changed 
significantly.  The  maximum  unit  flow  discharge  occurred  near  the  right 
bank  where  the  channel  was  the  deepest,  but  the  largest  flow  concentrations 
and  unit  bed-load  discharges  were  found  in  the  shallower  area  near  the 
left  bank. 

Figures  I-C-3.11,  I-C-3.12,  and  I-C-3.13  show  the  lateral  distri¬ 
butions  at  sections  11-1-1,  11-2-1,  and  12-1-1,  respectively.  The  width  of 
the  main  channel  decreases  considerably  between  sections  9-1  and  11-1, 
with  a  corresponding  increase  in  flow  depth.  The  figure  of  section  11-1-1 
sliows  that  depths  of  over  35  ft  were  measured  near  the  right  bank.  As 
seen  in  figure  I-A.l,  section  11-2  could  be  treated  as  three  smaller  sub¬ 
sections.  Flow  from  the  two  subsections  nearer  the  right  bank  reentered 
the  main  channel  upstream  from  section  12,  but  the  flow  through  the  sub- 
scx'.t ion  nearer  the  left  bank  continued  along  the  back  side  of  an  island 
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and  re-entered  the  main  channel  downstream  from  section  12.  The  largest 
depths,  mean  flow  velocities,  and  unit  water  and  suspended-sediment  dis¬ 
charges  at  section  11-2-1  were  measured  in  the  subsection  nearer  the  left 
bank,  but  the  largest  unit  bed-load  discharge  occurred  in  the  center  sub¬ 
section  and  thereafter  joined  the  main  channel  flow  upstream  from  section 
12-1.  Although  the  flow  quantities  varied  considerably  across  section 

11- 2-1,  the  median  bed-material  size  and  mean  concentration  were  practi¬ 
cally  constant.  The  width  of  the  main  channel  doubled  from  section  11-1 
to  section  12-1  which  was  the  widest  main-channel  section  of  the  study 
reach.  For  section  12-1-1,  shown  in  figure  I-C-3.13,  the  two  distinct 
peaks  in  the  unit  bed-load-discharge  distribution  each  represent  the  down¬ 
stream  effect  of  the  sediment  discharged  through  sections  11-1-1  and  11-2-1. 

Figures  I-C-3.14  through  I-C-3.26  for  sections  1-2-2  through 

12- 1-2,  respectively,  present  the  lateral  distributions  computed  from  the 
data  obtained  during  the  second  trip.  Water  discharges  for  the  second  trip 
were  low,  typically  one-third  of  those  during  the  first  trip.  Bed-load 
discharges  were  also  extremely  arall  and  therefore  large  percentage  dif¬ 
ferences  resulted,  which  were  subject  to  error. 

Figure  I-C-3.14  shows  the  lateral  distributions  at  section  1-2-2. 
Similar  to  the  first  trip,  the  largest  value  for  mean  suspended-sediment 
concentration  was  measured  near  the  right  bank.  This  concentration  was 
the  downstream  effect  from  the  confluence  of  the  DMR  and  MR.  The  mean 
sediment  concentration  in  the  DMR  on  the  previous  day  was  roughly  100  ppm. 
Another  possible  effect  on  the  sediment  concentration  and  bed-load  measure¬ 
ments  taken  near  the  right  bank  of  this  section  was  the  dredge-spoil- 
removal  operation  being  done  by  the  Corns  of  Engineers,  Rock  Island  District 
(COE(RI)).  Sand  was  being  taken  off  near  the  south  end  of  Fox  Island  by 
crane  and  transported  upstream  in  barges  to  L&D  19  for  cofferdam  construc¬ 
tion.  Minor  amounts  of  sand  being  spilled  into  the  river  may  have  affected 
the  measurements  taken  near  this  location. 

The  lateral  distributions  at  sections  2-1-2  and  3-2-2  are  shown 
in  figures  I-C-3.15  and  I-C-3.16,  respectively.  The  increases  in  the 
depth,  mean  velocity,  and  unit  bed-load  discharge  observed  near  the  left 
bank  of  section  2-1-2  compared  with  section  1-2-2  were  caused  by  the  flow 
between  the  small  downstream  island  ;md  Taylor  Island.  The  mean  concentra¬ 
tion  was  smaller  in  the  center  of  the  section,  as  was  ;»*en  lor  the  first 


q  xlO  q  *10 


3 


(ppm) 


q  xlO  q  xlO 


(ppm) 


XlO  q  XlO 


(uidd ) 


trip.  The  decreased  flow  velocities  in  the  shallow  center  of  the  section, 
which  reduced  the  suspended-sediment  transport  capacity,  and  the  downstream 
bifurcation  of  flow  at  section  3  caused  the  low  mean  concentration  observed. 
At  section  3-2-2,  the  larger  unit  bed- load  discharge  measured  near  the 
right  bank  was  produced  by  the  unstable  upstream  tip  of  Buzzard  Island. 

The  source  of  tlie  high  mean  concentration  observed  near  the  left  bank  was 
caused  by  a  flow  with  more  suspended  sediment  from  the  side  channel  behind 
Taylor  Island,  or  by  the  increased  sediment-transport  capacity  of  the  secon¬ 
dary  currents  formed  as  the  thalweg  crossed  the  river  between  sections  2-1 
and  3-2. 

Figures  I-C-3.17  and  I-C-3.18  present  the  lateral  distributions  at 
sections  4-1-2  and  5-1-2,  respectively.  The  highest  mean  suspended-sediment 
concentration  at  section  4-1-2  occurred  at  the  right  bank,  as  was  the 
case  in  the  first  trip,  due  to  the  combination  of  the  unstable  upstream 
tip  of  Buzzard  Island  and  the  downstream  influence  of  the  suspended-sediment 
inflow  from  the  MR.  At  section  5-1-2,  very  coarse  median  particle  sizes 
continued  to  be  found  roughly  400  feet  from  the  left  bank,  as  seen  in  figure 
I-C-3.18.  Close  to  the  left  bank  at  this  section,  where  the  flow  depth 
was  only  2.9  ft,  the  surface  flow  was  upstream  and  thus  gave  negative 
values  for  the  mean  flow  velocity  and  suspended-sediment  load. 

The  lateral  variations  at  sections  6-1-2  through  9-1-2  are 
presented  in  figures  I-C-3.19  through  I-C-3.22,  respectively.  This  area 
near  the  downstream  tip  of  Buzzard  Island  is  regarded  as  a  problem  shoaling 
area  that  requires  occasional  dredging.  The  figures  show  that  although 
the  magnitude  of  the  unit  flow  discharge  remained  nearly  constant  from  sec¬ 
tion  6-1-2  through  section  7-1-2,  the  maximum  unit  suspended-sediment  dis¬ 
charge  doubled,  while  the  unit  bed-load  discharge  remained  insignificant 
at  section  7-1-2.  The  increased  suspended-sediment  load  suggests  that 
scouring  was  occurring  between  these  two  sections.  Precise  interpreta¬ 
tion  is  difficult  since  the  measurements  at  these  two  stations  were  made 
on  different  days  (section  6-1-2  on  18  August  1978;  section  7-1-2  on  21 
August  1978).  Between  these  dates,  the  water  discharge  increased  by  an 
estimated  22  iiercont.  However,  rnc'asurements  at  sections  7-1-2  and  8-2-2 
wore  taken  on  the  same  day.  Ihe  figures  for  these  sections  reveal  (.hat  tin* 
maximum  unit  suspended-sediment  discharge  at  section  8-2-2  was  roughly 
one-half  of  its  maximum  value  at  section  7-1-2.  The  observed  reduction  in 
the  sediment-transport  capacity  of  the  river  is  attributable  to  the  stealing 
problem  observed  in  this  area. 
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There  are  three  general  features  of  the  MR  in  this  shoaling 
area  that  reduce  the  sediment-transport  capacity  of  the  flow.  First, 
the  river  width  increases  from  section  7-1  to  section  8-2,  which  increases 
the  cross-sectional  flow  area.  For  the  same  discharge,  the  increased  flow 
area  decreases  the  mean  flow  velocity,  resulting  in  a  decrease  in  the 
sediment-transport  capacity.  Consequently,  the  excess  sediment  is  deposited 
until  an  equilibrium  is  achieved  between  the  sediment-transport  capacity  and 
sediment  input.  Second,  the  tlialweg  crosses  the  river  from  the  left  bank 
at  section  5-1  to  the  right  bank  at  section  9-1.  Near  the  middle  of  this 
crossing  an  inflection,  or  so-called  ’'cross-over”,  results.  The  cross¬ 
over  region  generally  lacks  the  strong  secondary  currents  produced  by 
bends  in  the  river  that  significantly  increase  the  sediment -transport  capa¬ 
city  of  the  flow.  Finally,  the  bifurcation  of  the  flow  between  transects 
6  and  9  reduces  the  main  channel  flow  and  therefore  diminishes  the 
sediment  transport  capacity. 

Specific  features  of  sections  7-1-2  to  9-1-2  will  now  be  dis¬ 
cussed.  Figures  I-C-3.20  and  I-C-3.21  of  sections  7-1-2  and  8-2-2  show  that 
the  larger  depths  and  unit  water  discharges  occurred  near  the  left  bank 
of  Hunt  Island,  although  the  larger  mean  suspended-sediment  concentrations 
and  unit  suspended-sediment  discharges  occurred  over  the  navigation 
channel  near  the  right  bank  of  Buzzard  Island.  During  the  first  trip, 
the  larger  unit  suspended-sediment  discharges  were  measured  near  the 
left  bank  at  sections  7-1-1  and  8-2-1.  It  is  surmised  that  the  larger 
flow  discharges  during  the  first  trip  caused  a  greater  extent  of  lateral 
dispersion  of  suspended  sediment  around  Buzzard  Island  than  the  smaller 
discharges  during  the  second  trip.  The  two  peaks  seen  in  the  unit  bed-load - 
discharge  distribution  over  the  shallower  central  part  of  section  8-2-2 
were  due  to  the  effect  of  the  bifurcation  downstream  at  transect  9.  The 
figure  for  section  9-1-2  shows  that  the  largest  depths,  unit  water,  and 
suspended-sediment  discharges  occurred  along  the  navigation  channel  near 
the  right  bank.  Note,  though,  that  the  largest  mean  velocities  were  meas¬ 
ured  near  the  left  bank  where  the  flow  bifurcated  around  the  upstream 
tip  of  Huff  Island.  The  bed  load  samples  collected  at  this  section  were 
negligible.  Tabic  I-C-3.2  sunrnarizes  the  percentages  of  water  and  sediment 
discharges  that  flowed  into  the  side-channel  sections  7-2-2  and  9-2-2. 


42 


Section 


7-1-2 

7-2-2 


Table  I-C-3.2 

Water  and  Sediment  Discharges  at  Sections  7  and  9 


082178 

082178 


Q 

(cfs) 

44,509 

4.010 


(tons/ day) 

13,458 

467 


( tons/ day ) 


0.07 

%2/qbi=0^ 


9-1-2 

9-2-2 


082278 

082278 


30,703 

7,592 

Q2/Q1=25f; 


4,281 


0.3 

TO-®5 


The  lateral  distributions  at  sections  10-1-2,  11-1-2,  11-2-2, 
and  12-1-2  are  presented  in  figures  I-C-3.23  through  I-C-3.26,  respectively. 
The  figure  of  section  10-1-2  shows  that  the  largest  flow  quantities  were 
measured  near  the  right  bank  where  the  navigation  channel  was  located. 

At  section  11-1-2,  bed-material  samples  near  the  right  bank  consisted 
mainly  of  coarse  gravels,  while  samples  near  the  left  bank  were  mainly 
silt  and  clay.  The  figure  of  section  11-2-2  shows  that  the  largest  unit 
bed-load  discharge  (although  negligible)  was  measured  in  the  middle  sub¬ 
section  of  the  transect,  as  was  seen  for  the  first  trip.  However,  the 
mean  flow  velocities  were  slightly  greater  in  this  middle  subsection  instead 
of  in  the  subsection  nearest  the  left  bank,  as  was  the  case  for  the  first 
trip.  Note  tliat  the  figure  of  section  12-1-2  (figure  I-C-3.26)  shows 
that  the  bed-material  composition  near  the  right  bank  changed  appreciably 
from  the  first  trip  and  was  due  to  the  deposition  that  occurred  here  between 
trips,  which  will  be  discussed  later. 


4.  Longitudinal  distributions.  longitudinal  distributions  of 
the  major  flow  and  sediment  quantities  are  plotted  for  the  main-channel 
sections  in  figures  I-C-4.1  and  I-C.4.2,  and  for  the  complete  river 
cross  sections  in  figure  I -CM.  3.  The  longitudinal  variations  of  water 
(lisHniri?*,  ();  siis|xmdod-sediment  discharge,  Q  ;  and  Ixxi-load  discharge, 

Q  ;  arc  unvalued  in  figures  f-C-l.l  and  MV!. 3,  while  (he  variations  of 


(cfs)  (tons/ 


( tons/day )  ( tons/day ) 


Q  (xl'J 
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unit  water  discharge  q  (=Q/W),  mean  flow  velocity  U  (=<)/A),  unit  suspended- 

(^/W),  and  mean 
median  bed-mater ial  size,  are  presinted  in  figure  1-04.2. 

The  flow  and  sediment  discharges  for  the  conplete  river  width  were 
approximated  by  adding  the  flow  and  sediment  discharges  of  the  appropriate 
side-channel  sect  ion  (s)  to  the  discharges  measured  in  the  main-channel 
section.  Side  channels  were  chosen  according  to  their  location  and  the  day 
that  their  discharges  were  measured.  Since  the  flow  and  sediment  discharges 
estimated  for  the  complete  river  sections,  shown  in  figure  1-04.3,  vary 
little  from  the  discharges  in  the  main-channel  sections,  shown  in  figure 

1- 04 . 1 ,  the  following  discussions  will  focus  on  the  latter  figure,  with 
similar  results  implied  from  the  figure  I-C-^.3.  When  comparing  the  general 
trends  between  the  figures  1-04 . 1  and  1-04.3,  note  that  several  of  the 
scales  have  been  changed.  In  these  figures,  the  sections  are  equally  spaced 
on  the  horizontal  scale,  which  is  not  intended  to  represent  actual  distances 
between  stations.  Also,  sections  were  not  always  sampled  in  order:  samples 
from  section  11-1-1  were  collected  on  10  July  1978,  while  samples  frcm 
section  10-1-1  were  collected  on  11  July  1981.  Interpretation  of  the  longi¬ 
tudinal  variations  along  the  study  reach  is  difficult  since  the  MR  or  any 
other  natural  river  is  not  a  steady -state  system.  The  quantities  calculated 
and  plotted  for  each  section  were  measured  during  a  time  of  roughly  4  hours, 
while  the  complete  set  of  data  for  each  trip  was  collected  over  a  period  of 
approximately  2  weeks,  during  which  a  number  of  outside  influences  could 
alter  the  river  flow7  and  sediment  characteristics.  The  longitudinal  dis¬ 
tributions  will  be  discussed  next.  Note  the  use  of  subscripts  and  different 
scales  to  denote  the  two  trips  (e.g.  ,  and  Q^). 

For  the  first  trip,  a  sharp  increase  in  the  bed-load  discharge 
was  observed  between  sections  1-2-1  and  2-1-1.  The  measurements  were 
made  on  consecutive  days  (section  1-2-1  on  26  June  1978;  section  2-1-1  on 
27  June  1978).  Although  the  bed-load  discharge  increased,  the  water  and 
suspended-sediment  discharges  decreased.  As  seen  in  figure  I-A.l,  section 

2- 1  w ns  located  at  a  river  bend,  where  the  main  channel  width  was  narrower 
Ulan  at  sect  ion  1-2.  The  narrower  river  width  reduce*!  tlie  cress-section 
area  and  inenmsed  I  lie  mean  Now  velocity,  which,  l.o  a  minor  extent. 


sediment  discharge  qg  (=<3S/W),  unit  bed-load  discharge 
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increased  the  bed-load-transport  capacity.  The  strong  secondary  currents 
occurring  in  the  river  bend  are  the  primary  cause  of  the  bed-load  movement. 

The  entirely  different  cross-sectional  profiles  of  sections  1-2  and  2-1 
shown  in  figure  I-C-2.1  illustrate  the  effect  the  secondary  currents  have 
had  along  the  right  bank  of  section  2-1.  The  tributary  inflow  from  the 
Fox  River,  which  entered  the  MR  approximately  1.5  mi  upstream  from  section 
2-1,  might  have  also  contributed  to  the  large  bed-load  movement.  No  flow 
measurements  were  taken  in  the  Fox  River,  so  a  comparison  was  not  possible. 

The  bed-load  discharge  decreased  at  the  next  several  sections  downstream 
from  section  2-1-1. 

Heavy  rains  throughout  Iowa  during  the  first  trip  on  25  and  27 
June  1978  caused  the  increased  water  and  suspended-sediment  discharges  measured 
at  section  4-1-1  on  29  June  1978.  These  high  discharges  gradually  diminished 
over  the  following  week. 

Although  the  water  and  suspended-sediment  discharges  continued 

to  decrease  at  sections  6-1-1  and  7-1-1,  the  bed-load  discharge  displayed 

a  significant  rise  at  section  7-1-1  and  then  decreased  at  section  8-2-1. 

Table  I-C-4.1  surrmarizes  the  flow  and  sediment  discharges  measured  during 

the  first  trip  at  sections  6,  7,  and  8.  The  bed-load  discharge  increased  by 

a  factor  of  5  at  transects  6  and  7,  which  suggests  that  scouring  was  taking 

place  along  this  reach.  The  cross-sectional  profiles  of  section  7-1  in 

figure  I-C-2.1  indicate  that  some  scouring  had  occurred  here  between  trips. 

Although  roughly  16  percent  of  the  main  channel  flow  entered  the  side  channel 

(section  7-2)  and  the  river  width  increased  between  sections  6-1  and  7-1, 

o 

the  cross-sectional  area  decreased  from  an  estimated  39,200  ft  at  section 
2 

6-1-1  to  31,600  ft  at  section  7-1-1.  Estimation  of  the  percentage  of 
water  discharge  from  section  6-1  that  would  flow  through  section  7-1,  based 
on  the  first -trip  measurements,  gives  a  7  percent  increase  in  the  mean 
flow  velocity.  This  increased  mean  flow  velocity  would  increase  the  amount 
of  bed  load  being  transported  at  section  7-1,  although  not  to  the  extent 
that  was  observed  during  the  first  trip. 
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Table  I-C-4.1 

Water  and  Sediment  Discharges  at  Sections  G,  7,  and  8 


Section 

Date 

Q 

(cfs) 

(tons/day) 

% 

(tons/day) 

6-1-1 

070578 

123,949 

118,541 

414 

7-1-1 

070578 

101,830 

89,467 

1,958 

7-2-1 

070678 

16,155 

13,280 

102 

117,985 

102,747 

2,060 

8-2-1 

070678 

108,694 

75,436 

681 

Since  measurements  for  sections  7-1-1  and  8-2-1  were  taken  on  con¬ 
secutive  days,  flow  discharges  should  compare  favorably.  Although  the  flow 
discharge  increased  slightly  fran  7-1-1  to  8-2-1,  the  suspended-sediment  dis¬ 
charge  decreased  and  the  bed-load  discharge  at  8-2-1  was  only  one-third  of 
that  at  7-1-1.  The  mean  velocity,  U,  had  a  slight  reduction  frcm  3.22  ft/s 

n 

to  3.04  ft/s  as  the  cross-sectional  area  increased  from  31,600  ft*  to  35,700 
2 

ft  .  For  the  same  discharge,  the  mean  flow  velocity  would  decrease  roughly 
10  percent  from  sections  7-1  to  8-2.  This  reduced  mean  velocity  would  decrease 
the  sediment -transport  capacity  of  the  river,  thus  causing  the  deposition  of 
sediment  in  the  reach  between  sections  7-1  and  8-2.  This  area  located  near 
the  downstream  tip  of  Buzzard  Island,  as  discussed  earlier,  is  known  for  its 
frequent  shoaling  problems  (Naka/to  and  Kennedy ,  1977 ) . 

During  7-9  July  1978,  heavy  rain  and  flooding  in  northwest  Iowa 
caused  an  increase  in  flow  along  the  MR.  The  effect  of  this  rain  is  seen  in 
the  large  water  and  suspended-sediment  discharges  at  sections  10-1,  11-1,  and 
12-1.  Note  that  part  of  the  large  increase  observed  in  the  water  discharge 
between  sections  11-1  and  12-1  was  due  to  the  incoming  flow  through  section 
11-2.  Table  I-C-4.2  surrmarizes  the  water  and  sediment  discharges  through 
.sections  11  and  12  for  the  first  trip.  The  numbers  in  parentheses  indicate 
the  percentage  of  the  respective  total  discharge  passing  through  side  channel 
11-2. 

A  water  discharge  of  115,767  cfs  was  measured  on  11  July  197S 
at  section  10-1.  Since  the  water  discharge  measured  on  the  previous  clay 
at  sect  ion  11-1  was  approximately  2  jxavent  larger,  there  was  pract  ically 
iv*  change  in  I  lie  flow  condition::  1‘nni  It)  to  11  July  1978,  during  which 
transects  11  and  12  were  measured.  Then  I’ run  continuity  con: ;id< ‘rations, 
the  water  or  sediment  discharges  at  sect  ion  12-1-1  should  closely  approximate 
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the  discharge  measured  at  section  11-1-1  plus  part  of  the  incoming  dis¬ 
charge  from  section  11-2-1.  Approximately  54  percent  of  the  water,  sus¬ 
pended-sediment,  and  bed- load  discharges  measured  at  section  11-2-1  entered 
the  main  channel  upstream  from  section  12-1-1.  Using  this  54  percent  approxi¬ 
mation,  the  main-channel  water  discharge  at  section  12-1,  estimated  from 
section  11,  differed  by  less  than  4  percent  from  the  water  discharge  meas¬ 
ured  at  section  12-1.  However,  the  bed-load  discharge  measured  at  section 
12-1-1  was  less  than  that  measured  at  section  11-1-1,  without  considering 
the  extra  bed-load  discharge  passed  through  section  11-2-1  into  the  main 
channel.  In  addition,  an  estimated  49,800  tons/day  of  the  suspended-sediment 
discharge  at  section  11-2-1  entered  the  main  channel  upstream  from  section 
12-1  through  the  two  subsections  nearer  the  right  bank  of  section  11-2-1. 

When  added  to  the  suspended-sediment  discharge  through  section  11-1-1, 
the  total  exceeds  the  suspended-sediment  discharge  measured  at  section 
12-1-1  by  nearly  20  percent.  The  reduction  in  bed-load  discharge  resulted 
from  the  smaller  flow  velocity  at  section  12-1,  compared  to  those  at 
sections  11-1  and  11-2.  The  mean  velocity,  U,  shown  in  figure  I-C-4.2, 
decreased  from  3.71  ft/s  at  section  11-1-1  to  3.25  ft/s  at  section  12-1-1 
due  to  the  increased  river  width.  The  lateral  distributions  at  sections 
11-1-1,  11-2-1,  and  12-1-1  in  figures  I-C-3.11,  I -€-3.12,  and  I-C-3.13, 
respectively,  show  this  velocity  decrease  in  greater  detail  across  the 
transects.  The  lower  flow  velocities  across  section  12-1-1  diminished  the 
sediment-transport  capacity  of  the  river,  which,  therefore,  reduced  the 
suspended-sediment  and  bed-load  discharges.  The  cross-section  profiles 
of  section  12-1  in  figure  I-C-2.1  show  that  deposition  occurred  over  most 
of  the  led,  with  up  to  7-ft  rises  in  the  bed  elevation.  This  reach  of  the 
MR  is  sufficiently  wide  and  deep  that  large  amounts  of  sediment  deposition 
presently  do  not  critically  affect  the  navigation  channel.  According  to 
COE  (RI),  depths  of  over  60  ft  have  been  measured  in  the  main  channel  two 
miles  downstream,  at  RM  345. 

The  water  discharge  during  the  second  trip  sliowod  only  minor 
fluctuations,  as  scx'n  in  figures  I-C— 1.1  and  T-C— 1.2.  Note'  tliat  part  of 
the  drop  in  water  discharge  downstream  ftxmi  section  8-2-2  was  due*  to  the 
flow  bifurcation  at  transect  9.  Most,  of  the  I  low  through  the  sick'  ch;umo  1  , 
section  8-2-2,  reentered  the  main  channel  through  section  11-2-2  and 
caused  some  of  tlie  subsequent  rise  in  water  discharge  at  section  12-1-2. 
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Distinct  peaks  occurred  in  the  suspended-sediment  and  bed-load 
discharges  at  section  7-1-2.  This  coincidence  of  peaks  is  interesting, 
since  these  sediment  quantities  were  calculated  independently  of  each  other. 
This  peaking  occurred  in  the  same  shoaling  area  that  was  discussed  for 
the  first  trip  except  now  the  sediment  discharges  decreased  toward  trans¬ 
ect  9,  instead  of  transect  8.  Table  I-C^1.3  summrizes  the  measured 
flow  quantities  and  sampling  dates  for  sections  6  through  9. 


Table  I-C-4.3 


Water  and  Sediment  Discharges  at  Sections  6  through  9 


Section 

Date 

Q 

(cfs) 

Qs 

(tons/day) 

% 

(tons/day) 

6-1-2 

081878 

39,790 

4,046 

3.2 

7-1-2 

082178 

44,509 

13,458 

30.0 

7-2-2 

082178 

4,010 

467 

0.07 

48,519 

13,925 

30.07 

8-1-2 

082278 

1,363 

243 

0.3 

8-2-2 

082178 

46,400 

8,523 

21.0 

47,763" 

8,766 

21.3 

9-1-2 

082278 

30,703 

4,281 

1.4 

9-2-2 

082278 

7,592 

669 

0.3 

38,295 

4,950 

1.7 

Due  to  the  changes  in  flow  conditions  between  the  days  when 
sections  6  and  7  were  measured,  the  primary  cause  of  the  increased  sedi¬ 
ment  discharge  observed  at  section  7  cannot  be  resolved.  Since  measure¬ 
ments  at  sections  7-1-2  and  8-2-2  were  completed  on  the  same  day,  the  water 
and  sediment  discharges  should  compare  closely.  The  4  percent  difference 
seen  in  the  water  discharge  was  due  to  inaccuracies  of  the  measurements 
and  any  minor  flow  fluctuations  that  occurred  during  the  8-hr  sampling 
period.  The  significant  reduction  in  the  suspended-sediment  discharge 
illustrates  the  decreased  sediment-transport  capacity  of  the  flow  at  sec¬ 
tion  8-2-2.  The  overnight  decrease  in  flow  conditions  from  21  to  22  August 
1978  is  believed  to  have  caused  the  diminished  sediment  discharges  that 
were  measured  at  transect  9. 
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The  longitudinal  profile  of  the  composite  median  diameter,  £Lq, 
shows  that  the  particle  sizes  were  generally  finer  during  the  second  trip. 
This  decrease  in  was  likely  due  to  the  frequent  rains  throughout  the 
Midwest  during  June  and  July,  which  produced  high  water  discharges  and  flow 
velocities  in  the  MR  that  transported  large  amounts  of  fine  sediment.  These 
high  flow  conditions  were  evident  during  the  first  trip.  The  second 
trip  in  August,  however,  found  water  discharges  and  flow  velocities  roughly 
one-third  and  one-half ,  respectively,  of  those  measured  during  the  first 
trip.  This  reduction  in  flow  significantly  decreased  the  river’s  sediment- 
transport  capacity  so  that  the  suspended-sediment  discharges  during  the 
second  trip  were  only  approximately  4  percent  of  those  during  the  first 
trip.  Consequently,  part  of  the  sediment  that  was  in  suspension  during  the 
first  trip  settled  out  of  the  flow  and  became  bed  material.  The  sediment 
deposited  from  the  lower  flows  was  finer  than  the  median  particle  size  in 
the  existing  bed,  which  resulted  in  a  higher  proportion  of  fine  sediment, 
and  therefore  a  finer  median  bed-material  size. 

One  final  observation  was  that  the  composite  median  diameter,  D^, 
at  section  12-1  was  larger  during  the  second  trip  than  the  first  trip. 

This  observation  is  in  agreement  with  the  previous  discussion,  since  the 
bed  load  deposited  at  section  12-1  came  from  upstream  reaches  where  the 
mean  particle  size  was  larger. 

5.  Velocity  and  suspended-sediment-concentration  profiles. 
Velocity  and  concentration  profiles  were  constructed  for  the  DMR  and  the 
main-channel  sections  of  the  MR  where  detailed  flow  velocity  and  suspended- 
sediment  -concent  rat  ion  measurements  were  taken.  The  method  of  least 
squares  was  used  to  calculate  the  equation  of  the  line  representing  each 
velocity  and  suspended-sediment-concentration  profile.  Representative 
profiles  are  shown  in  figure  l-C-5.1  for  the  vertical  141  m  from  the  right 
bank  of  section  9-1-1. 

The  velocity  profiles  were  represented  by  the  well-known  loga¬ 
rithmic  relation: 


n 


log  (y/d )  i  <  < >n: ;l  . 


(1.1  ) 


whore  u  *  =  /gdS ; 


is  the  gravi lul  ional  onusUirt.;  S  is  the  energy  sin]**; 


and k  is  the  Karman  constant.  The  quantity  2.3  u*/c  was  determined  for 
each  velocity  profile  as  the  increment  in  velocity  over  a  logarithmic 
cycle  of  y.  The  shear  velocity,  u+,  was  computed  using  the  water-surface 
slope  between  the  upstream  and  downstream  pool  elevations  at  Keokuk,  Iowa, 
and  Canton,  Missouri,  respectively,  and  the  depth  of  flow  at  the  vertical. 
Finally,  k  was  obtained  from  the  calculated  values  of  2.3u+/k  and  u*. 

The  velocity  profiles  of  the  first  trip  exhibited  a  steeper 
slope  than  the  velocity  profiles  of  the  second  trip.  The  energy  gradients 
of  the  first  trip  were  greater  than  those  of  the  second  trip  by  a  factor 
of  approximately  5.  Since  the  mean  flow  depth  in  the  pool  was  maintained 
fairly  constant  by  the  downstream  gate  operation,  the  larger  energy  gra¬ 
dients  caused  larger  values  of  u+.  Consequently,  the  quantity  2.3u*/<  in 
(1.1)  was  larger  which  caused  the  velocity  profiles  of  the  first  trip  to 
be  inclined  more  than  those  for  the  second  trip. 

The  concentration  profiles  were  represented  by  the  Rouse  equa¬ 
tion: 


C_  _  rd-y  a  .z 
r  y  v  d-a' 


where  z  = 


"a 

w 


:>KUa 


(1.2) 

is  the  Rouse  number;  w  is  the  particle  fall  velocity;  $  is 


the  ratio  between  the  sediment  and  momentum  turbulence  exchange  coeffi¬ 
cients;  and  represents  the  suspended-sediment  concentration  at  9ome 

a 

reference  level,  a,  above  the  river  bed.  The  values  of  the  Rouse  number 
determined  from  the  concentration  profiles  of  the  first  and  second  trips 
are  tabulated  in  table  I-C-4.4.  As  can  be  seen  in  the  table,  the  values 
of  z  were  generally  smaller  during  the  first  trip  than  the  second  trip. 

The  smaller  values  of  z  that  resulted  from  the  near-vertical  concentration 
profiles  indicate  that  the  suspended-sediment  concentration  was  uniformly 
distributed  over  the  deeper,  highly  turbulent  flows  of  the  first  trip. 

All  but  two  of  the  concentration  profiles  had  positive  values  for  z,  which 
indicate  the  generally  observed  increase  of  suspended-sediment  concentration 
from  the  water  surface  to  the  river  bed  (increasing  (d-y)/y).  The  two 
nognt  ive  values  calculated  for  z  resulted  from  the'  least-squares  analysis 
of  inexact  susjx.'ndcxl-sediment  concentration  measurements,  which  were  less 
than  50  ppm  at  both  sections  and  varied  less  than  7  ppm  over  the  flow 
depth.  The  table  also  lists  the  values  of  the  Karman  constant,  <  , 


Table  I-C-4.4 


Calculated 

Values  of  u^,  <  ,  C,  and  z  for  the  Major  Verticals 

Section  and 

Trip  No. 

Distance  from  k  C  z 

R.B.  (m)  (ft/s)  (ppm) 

13-1-1 

99 

0.241 

0.30 

1791 

0.0076 

1-2-1 

162 

0.170 

0.26 

621 

0.032 

2-1-1 

174 

0.168 

0.31 

218 

0.037 

3-2-1 

353 

0.187 

0.42 

142 

0.161 

4-1-1 

183 

0.211 

0.44 

1193 

0.018 

5-1-1 

260 

0.210 

0.52 

690 

0.015 

6-1-1 

206 

0.206 

0.26 

428 

0.016 

7-1-1 

111 

0.131 

0.46 

431 

0.0081 

8-2-1 

151 

0.179 

0.53 

258 

0.051 

9-1-1 

141 

0.186 

0.55 

292 

0.030 

10-1-1 

36 

0.191 

0.32 

422 

0.015 

11-1-1 

30 

0.230 

0.29 

679 

0.0094 

12-1-1 

294 

0.199 

0.40 

465 

0.0018 

13-1-1 

80 

0.235 

0.98 

358 

0.024 

13-1-2 

86 

0.145 

0.82 

101 

0.175 

1-2-2 

170 

0.075 

0.25 

35 

0.036 

2-1-2 

56 

0.072 

0.32 

84 

0.163 

3-2-2 

262 

0.089 

0.21 

21 

0.126 

4-1-2 

305 

0.073 

0.46 

26 

0.081 

5-1-2 

246 

0.078 

0.66 

26 

0.063 

6-1-2 

188 

0.078 

0.60 

31 

0.038 

7-1-2 

119 

0.080 

0.40 

182 

0.013 

8-2-2 

141 

0.075 

0.46 

147 

0.073 

9-1-2 

124 

0.077 

0.35 

72 

0.041 

10-1-2 

44 

0.078 

0.25 

69 

0.168 

11-1-2 

19 

0.092 

0.65 

31 

-0.012 

12-1-2 

291 

0.075 

0.35 

47 

-0.0026 
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calculated  from  the  logarithmic  velocity  profiles,  and  the  mean  suspended- 
sediment  concentration,  C,  determined  from  depth- Integrated  suspended- 
sediment  samples.  These  values  exhibit  a  large  amount  of  scatter,  espe¬ 
cially  for  the  analler  concentrations  observed  during  the  second  trip  and 
for  those  points  representing  the  DMR.  Computed  values  of  <  ranged  from 
0.21  to  0.98,  while  for  clear  fluids,  the  value  of  <  is  approximately  0.4. 

One  last  observation  was  that  the  calculated  values  of  z  at 
section  7-1  for  both  trips  were  consistently  lower  than  those  at  nearby 
sections.  Lower  values  of  z  signify  more  uniformly  distributed  suspended- 
sediment-concentration  profiles.  Note  that  section  7-1  also  had  consis¬ 
tently  high  bed-load  discharges,  as  seen  in  figures  I-C-4.1  and  I-C-^4.2. 

This  section  is  just  upstream  from  the  problem  shoaling  region  near  the 
downstream  tip  of  Buzzard  Island. 

6.  Flow  and  sediment -discharge  relationships.  Figures  I-C-6, 1 
and  I-C-6. 2  illustrate  the  empirical  power-law  relationships  that  were 
formulated  between  the  unit  bed-load  discharge,  q^,  and  water  discharge, 

Q,  using  the  method  of  least  squares.  Figure  I-C-6. 1  presents  the  correla¬ 
tion  between  q~  and  U  for  the  main-and  side-channel  sections  along  the  .MR 
study  reach.  Note  that  only  sections  with  values  of  qQ  greater  than  10 
tons/ft/day  were  included.  Although  the  side-channel  sections  were  gen¬ 
erally  much  smaller  in  cross-sectional  area  and  had  smaller  water  and 
sediment  discharges  than  the  main-channel  sections,  the  unit  bed-load  dis¬ 
charges  and  mean  flow  velocities  for  both  the  main-  and  side-channel  sec¬ 
tions  were  of  comparable  magnitude.  The  result  for  the  least -squares  analy¬ 
sis  yielded 

q_  =  6.36  x  10"4  U5'5  (1.3) 

D 

Figure  I-C-6. 2  illustrates  the  relationship  between  q^  and  U 
for  the  main-channel  sections  (including  section  11-2-1),  which  is  repre¬ 
sented  by 


q  =  5.66  x  10-4  U5'7 


(1.4) 


This  ivuul  h>n  shews  1 1 1;  1 1  q  is  v<  *r\  sensitive  1«>  'imnirs  in  I  ,  ;«nd  licit  ;i 
peiveul  inrl*e;ise  ill  lilt1  !  I  n\\  \  <  •  h  »e  i  I  y  W'llld  lx-  < *\|  x  ■<  *1  » *<  I  l«»  producf 
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approximately  a  C -percent  increase  in  the  unit  bed-load  discharge.  The 
expected  increase  in  the  unit  bed-load  discharge  for  an  increase  in  the 

mean  flow  velocity  can  be  estimated  from  (1.4).  For  example,  closing 
side  channel  9-2  would  increase  the  mean  flow  velocity  through  section 
9-1  by  an  estimated  28  percent  (the  average  of  the  values  of  31  percent 
during  the  first  trip  and  25  percent  during  the  second  trip,  as  listed 
in  tables  I-C-3.1  and  I-C-3.2,  respectively),  assuming  the  cross-sectional 
area  remained  constant.  From  (1.4),  this  increased  flow  velocity  would 
be  expected  to  increase  the  unit  bed- load  discharge  in  section  9-1  by 
a  factor  of  approximately  2.6.  Similarly,  the  closure  of  section  7-2 
would  increase  U  in  section  7-1  by  an  estimated  12  percent  (the  average 
of  the  values  of  16  percent  during  the  first  trip  and  9  percent  during 
the  second  trip,  as  shown  in  tables  I-C-3.1  and  I-C-3.2,  respectively), 
assuming  a  constant  cross-sectional  area.  This  increased  U  would  be 
expected  to  increase  the  bed-load  discharge  through  section  7-1  by  approxi¬ 
mately  68  percent. 

Figure  I-C-6.3  shows  the  correlation  between  the  suspended-sedi¬ 
ment  discharge,  Q  ,  and  the  water  discharge,  Q,  for  the  MR  main-channel 
s 

sections  (note  that  section  11-2  was  included  due  to  its  large  flow  dis¬ 
charge).  The  power-law  relation  is 

Q  =  1.44  x  lCf8  Q2,5  (1.5) 

s 

The  suspended-sediment  discharge  is  not  as  sensitive  to  clianges  in  the 
water  discharge  as  the  unit  bed-load  discharge  is  to  clianges  in  the  mean 
flow  velocity.  It  is  shown  by  (1.5)  that  for  the  28  percent  average 
increase  in  Q  expected  if  section  9-2  were  closed,  the  suspended-sediment 
discharge  through  section  9-1  would  increase  by  70  percent.  The  12  per¬ 
cent  average  increase  in  Q  through  section  7-1  expected  if  section  7-2 
were  closed  would  increase  the  suspended-sediment  discharge  through  section 
7-1  by  approximately  30  percent.  Therefore,  the  closure  of  section  7-2 
or  9-2  would  signi  f  icantly  increase  both  the'  bed-load  and  susjx'nded- 
.stxliiiionf-  discharges  thiongh  the  main  eluinnel .  According! y ,  it  can  lx* 
concluded  that  the  closure  ot  either  side  send  ion  would  help  alleviate 
the  slioaling  problem  near  the  downstream  tip  ol  Buzzard  Island. 


Q  (tons/doy) 


F  iy»urc! 


hU 


Q (cfs) 


T-C-6.3  Relationship  between  0  and  O  for  the  main-channel 

s 


sections 
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One  final  point  worth  mentioning  is  the  portion  that  the  bed¬ 
load  discharge  contributed  to  the  measured  total -load  discharge,  = 

Qg  +  Qg.  From  the  data  collected  during  this  study,  it  appeared  that  in 
all  but  a  few  sections  the  bed-load  discharge  constituted  less  than  1  per¬ 
cent  of  the  total  measured  sediment  discharge,  especially  for  the  low  dis¬ 
charges  during  the  second  trip.  Note  that  the  suspended-sediment  load, 
was  considered  to  be  primarily  wash  load. 

D.  Conclusions.  The  principal  conclusions  derived  from  the 
present  study  are  sunmarized  as  follows: 

1.  Ibere  are  three  general  features  of  the  MR  near  the  shoaling 
reach  at  Buzzard  Island  to  which  the  shoaling  problem  is  attributed.  First, 
the  main- channel  flow  bifurcates  at  two  locations  in  this  reach.  Second, 
the  river  widens  and  reduces  the  sediment -transport  capacity  since  the 
increased  cross-sectional  flow  area  decreases  the  mean  flow  velocity.  Finally, 
located  near  the  downstream  tip  of  Buzzard  Island  is  a  cross-over  reach 
which  generally  lacks  the  strong  secondary  currents  that  are  produced  by 
bends  in  the  river  and  significantly  increase  the  sediment -transport  capa¬ 
city  of  the  flow  there. 

2.  The  sediment  causing  the  shoaling  problem  originates  primarily 
from  the  DMR.  The  smaller  energy  slope  and  velocity  of  the  MR,  in  conpari- 
son  with  those  of  the  DMR,  are  unable  to  transport  the  coarse  sediment 
inflow  from  the  DMR,  and  deposition  results. 

3.  During  the  high  discharges  typical  of  the  first  trip,  the 
large  flow  velocities  in  the  MR  swept  the  sediment  inflow  from  the  DMR 
abruptly  downstream  from  the  confluence.  Since  the  river  exhibits  little 
curvature  downstream  from  the  confluence,  the  cross-stream  mixing  was  grad¬ 
ual.  Consequently,  the  mean  suspended-sediment  concentrations  measured 
at  the  downstream  sections  were  considerably  higher  near  the  right  bank, 
and  continued  so  for  roughly  12  mi.  In  addition,  the  first-trip  measure¬ 
ments  reveal ed  that  the  right  bank  suspended-sediment  concentrations  0  to 
10  mi  downstream  from  the  confine  'nee  were  roughly  one-third  of  the  concen¬ 
trations  measured  in  the  DMR. 


1  j 
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4.  The  velocity  profiles  of  the  major  verticals  along  the 
study  reach  were  described  well  by  the  logarithmic  relation,  given  by 
(1.1).  The  Karman  constant,  k  ,  was  found  to  vary  widely  for  concentra¬ 
tions  less  than  1,800  ppm,  when  plotted  against  the  mean  suspended- 
sediment  concentrations  measured  at  the  respective  verticals.  This 
variation  was  especially  true  for  the  DMR  and  the  small  concentrations 
found  in  the  MR  during  the  second  trip.  The  suspended-sediment -concen¬ 
tration  profiles  for  the  major  verticals  were  generally  well-described  by 
the  Rouse  distribution,  given  by  (1.2). 

5.  Measurements  revealed  a  significant  increase  in  the  bed¬ 
load  discharge  at  section  7-1,  upstream  from  the  shoaling  area,  with  a 
subsequent  decrease  in  the  bed-load  discharge  at  section  8-2  immediately 
downstream  from  the  shoaling  area.  It  was  also  found  that  the  river  bed 
had  been  scoured  at  section  7-1  between  the  time  that  the  twa  trips 
were  conducted. 

6.  Power-law  relations  formulated  for  the  main-channel  sections, 
given  by  (1.4)  and  (1.5),  show  that  the  unit  bed-load  discharge,  qn,  varies 

O 

as  the  5.7-power  of  the  mean  flow  velocity,  U,  while  the  suspended-sediment 
discharge,  Qg,  varies  as  the  2.5-power  of  the  water  discharge,  Q. 

7.  Flow  measurements  showed  that  approximately  12  percent  of 
the  main-channel  flow  entered  the  side  channel,  section  7-2,  while  more 
than  25  percent  of  the  main-channel  flow  entered  the  side  channel,  section 
9-2.  This  flow  bifurcation  reduces  the  main-channel  mean  velocity,  which 
in  turn  diminishes  the  sediment -transport  capacity  and  causes  sediment 
deposition.  Application  of  the  power-law  relation  leads  to  estimates  that 
closure  of  section  7-2  would  increase  the  bed-load  discharge  and  suspended- 
sediment  discharge  through  section  7-1  by  approximately  70  percent  and 

30  percent,  respectively.  Similarly,  the  closure  of  section  9-2  would  in¬ 
crease  the  bed- load  discharge  by  a  factor  of  2.6,  and  increase  the  suspended- 
sediment  discharge  by  roughly  70  percent.  Hence,  it  would  appear  tliat 
closure  of  either  section  7-2  or  9-2  would  significantly  help  in  alleviating 
the  si  km.  I  mg  problem.  Note  that  the  estimates  Tor  IkH.Ii  .‘iortioiis  were  ob  ¬ 
tained  by  assuming  that  the  cross - :  :ec|  iotial  flow  area:;  lor*  seel  ions  V  I  and 
P-2  would  matin  relal  ivrly  iinisla.nl  a  Her  closure. 
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8.  The  effect  that  the  increased  sediment  discharge  would 
have  on  the  reach  downstream  from  the  shoaling  area  after  the  closure 
of  the  side  channels  appears  minor.  In  spite  of  the  low  flow  regime  during 
the  second  trip,  depths  measured  along  the  navigation  channel  at  section 
12-1  were  well  over  25  ft,  and  thus  posed  no  problem  to  navigation.  More¬ 
over,  flow  depths  of  over  60  ft  have  been  recorded  in  the  main  channel 
by  00E  (RI)  2  mi  farther  downstream,  at  HM  345.  Thus,  it  would  appear 
that  the  extra  sediment  being  transported  and  possibly  deposited  through 
this  downstream  reach  as  a  result  of  the  closure  of  sections  7-2  or  9-2 
would  not  critically  affect  the  navigation  channel. 
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II.  AN  EVALUATION  OF  POUR  NUMERICAL  MODELS 


Pool  20  of  the  MR,  extending  upstream  from  L&D  20  at  Canton, 
Missouri,  to  L&D  19  at  Keokuk,  Iowa,  was  the  site  of  two  recent  field 
studies  conducted  by  the  IIHR  (Nakato  and  Kennedy,  1977;  Vadnal,  1979). 
Localized  areas  of  shoaling  near  Fox  Island,  RM  355-56,  and  Buzzard  Is¬ 
land,  RM  349-50,  presented  troublesome  conditions  for  barge  navigation  due 
to  the  shallow  water  depths  that  resulted.  A  mathematical,  computer-based 
model  able  to  simulate  the  flow  conditions  in  this  reach  would  be  extremely 
useful  for  the  prediction  of  future  deposition  and  scour  trends,  as  well 
as  for  a  comparison  of  the  predicted  consequences  resulting  from  the  imple¬ 
mentation  of  various  proposed  corrective  measures.  Numerous  mathematical 
models  representing  the  state  of  the  art  are  currently  available.  This 
section  of  the  report  evaluates  the  performance  of  four  numerical  models 
that  were  applied  to  the  Pool  20  reach  of  the  MR  (RM  343.2  -  364.2).  These 
models  include:  the  HEC-6  model  developed  by  the  Hydrologic  Engineering 
Center,  U.S.  Army  Corps  of  Engineers;  the  UUWSR  and  SUSR  models  developed 
at  Colorado  State  University;  and  the  CHAR2  model  developed  by  Sogreah,  a 
consulting  firm  in  Grenoble,  France. 

A.  The  HEC-6  Mathematical  Model. 

1.  General  model  description.  The  HEC-6  mathematical  model  was 
supplied  by  the  Hydrologic  Engineering  Center  (HEC),  U.S.  Army  Corps  of 
Engineers.  The  model  is  a  one-dimensional  steady-flow  simulation  program 
designed  to  analyze  scour  and  deposition  in  rivers  and  reservoirs.  Cross 
sections  are  each  subdivided  into  a  part  which  has  a  movable  bed,  and  a 
Part  with  a  fixed  bed.  The  entire  movable-bed  portion  moves  vertically, 
due  to  degradation  and  aggradation,  while  the  other  part  remains  fixed 
throughout  the  simulation.  The  model  cannot  simulate  the  development  of 
m  nders,  the  lateral  distribution  of  sediment  load  across  a  cross  section, 
o.  density  and  secondary  currents.  Bed  forms  are  considered  only  indirectly 
by  varying  Manning's  roughness  coefficient  (n)  with  water  disci  large  or 
stage  elevation.  Additional  Features  of  the'  model  include  its  ability  to 
account  for  sediment  part  icle  armoring  and  dredging  ojX'rat.  ions. 


69 


The  Initial  IIHR  field  study  took  place  from  I»Jay  through  September 
1976,  while  the  corrplementary  study  occurred  from  June  through  August  1978. 
Consequently,  the  data  from  the  1973  study  provided  the  basis  for  evaluating 
the  results  sought  from  a  mathematical  program  that  utilized  the  1976  data 
for  its  initial  conditions.  Thus,  the  model-study  procedure  first  required 
the  conversion  of  the  information  gathered  from  the  1976  field  study  into 
data  that  were  us  id  for  the  initial  conditions  of  the  model.  The  HEC-6 
simulation  model  program  was  then  run  using  monthly-,  weekly-,  and  daily- 
averaged  flow  quantities  for  the  28-month  period  from  May  1976  through 
August  1978.  The  results  of  the  three  runs  are  discussed  and,  where  pos¬ 
sible,  compared  to  the  information  gathered  during  the  1978  field  study. 
Finally,  an  evaluation  of  the  applicability  of  the  HEC-6  mathematical  model 
for  predicting  deposition  and  scour  trends  is  made. 

Figure  II-A-1.1  shows  a  detailed  index  map  of  the  study  reach  along 
the  MR.  Figure  II-A-1.2  presents  a  schematic  model  outline  of  the  reach. 
Twenty-seven  cross  sections  were  included  over  the  21-mile  reach  between 
the  upstream  boundary  at  Keokuk,  Iowa  (LSD  19)  and  the  downstream  boundary 
at  Canton,  Missouri  (L&D  20).  Note  that  the  tributary  entry  point  for 
the  DMR  is  at  RM  361.4. 

The  description  of  the  data  required  for  the  model  is  presented 
in  section  2,  while  section  3  describes  the  model  calibration.  Section  4 
explains  several  sensitivity  tests  that  were  made,  and  section  5  presents  the 
results  of  the  simulation  runs. 

2.  Data  availability. 

a.  Geometric  data.  Cross  sections  within  the  problem-shoaling 
areas  near  Fox  and  Buzzard  Islands  were  obtained  from  the  1976  field  study. 
Other  cross  sections  were  obtained  from  topographic  maps  of  the  area  that 
were  compiled  in  a  1945  survey  by  the  OOE(RI).  The  cross  sections  taken 
directly  from  the  1976  field-study  report  were  more  detailed  and  included 
the  cross-sectional  areas  of  adjacent  side  ctumnols.  Reach  lengths  between 
cross  sections  vary  from  0.16  mi  in  one  of  the'  shoaling  areas  to  1.2  mi  at 
the  upstream  mxlel  boundary,  with  a  t  ypical  reach  length  of  1.0  ml .  The 
movab Le-bed  [x>rt ion  of  each  of  the  field-study  cross  sections  was  chosen 
to  correspond  to  that  portion  of  the?  bed  with  a  bed-material  particle- 
size  distr ibution  primarily  in  the  sand  rangf'.  For  other  cross  sections, 
the  rrovab !  o-txnr  J  fx>rtion  was  chosen  on  I  ho  basis  of  engineering  .judgment  and 
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b.  Sediment  data.  The  sediment  data  required  for  the  model  con¬ 
sist  of  the  boundary  inflowing  sediment  loads  and  their  size  gradations 
for  the  DMR  and  MR,  along  with  the  bed-material  size  gradation  for  each 
cross  section.  A  10-year  record  of  water  and  suspended-sediment  discharges 
measured  directly  downstream  from  L&D  19  was  used  to  formulate  a  logarithmic 
water-sediment- in  flow  relationship.  Figure  II-A-2.1  shows  the  least- 
squares  fit  of  the  monthly-averaged  data  and  the  derived  transport  rela¬ 
tionship: 


_o  p  44 

0^  =  2.39  xlO°  Q  (2.1) 

in  which  Qg  (tons/day)  is  suspended- load  discharge  and  Q(cfs)  is  water 
discharge.  Table  II-A-2.1  presents  the  breakdown  of  the  suspended-sediment 
load  into  standard  size  classifications  and  percentages  of  total  load  for 
the  MR.  Note  that  only  one  size  fraction  of  silt  was  considered,  although 
four  sizes  are  available  in  the  HEC-6  program.  Transport  relationships 
were  then  formulated  for  each  size  fraction.  Since  the  upstream  boundary 
of  the  model  is  directly  downstream  from  a  lock  and  dam,  the  bed- load, 
discharge  was  assumed  to  be  negligible  due  to  the  settling  effect  of  par¬ 
ticles  in  the  backwater  area  approaching  the  dam.  In  general,  the  sus¬ 
pended  load  of  the  MR  in  the  study  reach  consists  almost  entirely  of  wash 
load.  The  sediment-discharge  records  used  are  considered  to  be  reliable, 
and  thus,  the  equation  from  the  least -squares  fit  was  applied  throughout 
the  study  without  modification. 

Table  II-A-2.1 

Particle-Size  Distribution  as  Percent  of  Suspended  Load 
(Mississippi  River  at  Keokuk,  IA) 

Classification  and  Size  Percent  of  Suspended  Load 

Clay 

(<  0.004  rrm)  25.6 

Silt  (medium) 

(0.016-0.031  rrm)  73.8 

Very  Fine  Sand 
(0.0625-0.125  rrm) 


0.0 


l-’i^uru  1  1  -A-J  .  1  Relationship  between  ()  aiul  (}  for  the  Mississippi  River 
at  Keokuk,  Iowa 
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Sediment -discharge  records  for  the  DMR  near  Keokuk,  Iowa,  are 
not  as  complete  as  those  records  for  the  MR.  Therefore,  only  an  approxi¬ 
mate  water-sediment -discharge  relationship  was  formulated  from  four  months 
of  daily  suspended-sediment  measurements  taken  at  St.  Francisville ,  Mis¬ 
souri,  by  the  USGS.  The  initial  water-sediment-discharge  relation  obtained 
from  a  least-squares  analysis  is 

Qs  =  5.284  x  10-5  Q2,10  (2.2) 

Similarly,  a  least-squares  analysis  of  the  bed-load  discharge  versus  the 
water  discharge  gives  the  relation: 

Qg  =  1.19  x  10-8  Q2,59  (2.3) 

in  which  (tons/day)  is  the  bed-load  discharge.  Table  II-A.2.2  shows 
the  size  fractions  and  their  percentages  of  the  suspended-sediment  dis¬ 
charge  and  bed-load  discharge  for  the  DMR,  which  are  of  considerable  Impor¬ 
tance  when  analyzing  the  sediment-transport  capacity  of  the  i’IJR  downstream 
from  the  confluence  of  the  two  rivers.  The  larger  energy  gradient  of  the 
DMR  permits  it  to  transport  much  coarser  particles  than  the  MR,  as  can  be 
seen  by  comparing  tables  II-A-2.1  and  2.  The  larger-sized  particles 
that  are  discharged  into  the  MR  from  the  DMR  are  consequently  deposited  in 
the  areas  where  other  features  of  the  MR,  such  as  cross-overs,  interact  to 
decrease  the  sediment-transport  capacity  of  the  flow. 

As  mentioned  previously,  the  sediment  data  on  which  the  two 
DMR  relationships  for  the  suspended-  and  bed-load  discharges  are  based 
do  not  represent  corrplete  records,  and  so  the  derived  relationships  are 
questionable  and  may  require  modification  to  better  simulate  conditions 
observed  during  the  field  studies. 

The  bed-material  composition  of  each  cross  section  was  obtained 
either  directly  from  the  1976  field  study,  or  by  assuming  that  nearby 
cross  sections  had  identical  compositions  to  those  cross  sections  surveyed 
in  the  1976  study.  The  five  sections  from  RM  343.2  to  RM  317.0  used  the 
bed-material  compositions  from  the  1978  study,  since  the  1976)  .study  reach 
did  not  extend  that  far  downstre;un. 
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Table  I  l-A-2.2 

Particle-Size  Distribution  as  Percent  of  Suspended  and  Bed  Loads 

(Des  Moines  River) 


Classification  and  Size 

Percent  of  Suspended  Load 

C  lay 

15.0 

(<  0.004  mm) 

Silt  (med i urn) 

60.0 

(0.016-0.031  mm) 

Very  Fine  Sand 

25.0 

(0.0625-0. 125  mm) 

Percent  of  Bed  Load 

Fine  Sand 

1  .6 

(0. 125-0.250  mm) 

Med i urn  Sand 

8.0 

(0.25-0.50  mm) 

Coarse  Sand 

50.0 

(0.5-1 .0  mm) 

Very  Coarse  Sand 

25.6 

( 1 .0-2 .0  mm) 

Very  Fine  Gravel 

10.2 

(2. 0-4.0  mm) 

rine  Gravel 

4.6 

(4. 0-8.0  mm) 

76 


c.  Hydrologic  data.  The  required  hydrologic  data  consist  of 
the  MR  water-discharge  hydrograph,  the  water-surface  elevation,  and  the 
water  temperature  at  L&D  20.  This  location  and  its  flow  conditions  repre¬ 
sent  the  downstream  boundary  and  the  downstream  boundary  conditions  of  the 
model.  In  addition,  the  water-discharge  hydrograph  of  the  MIR  must  be 
known.  This  hydrograph  was  used  as  a  boundary  condition  for  the  tributary 
inflow  point. 

The  water-discharge  hydrographs  of  the  MR  at  L&D  20  and  the 
DMR  were  employed  to  compile  the  monthly-,  weekly-,  and  daily-averaged 
flow  conditions  that  were  used  for  the  model  study.  The  same  averaging 
was  done  for  the  water-surface  elevations  at  L&D  20.  However,  only  monthly- 
averaged  water  temperatures  were  employed  during  the  study.  The  tempera¬ 
tures  are  necessary  for  the  computation  of  sediment  particle  fall  veloc¬ 
ities. 

3.  Model  calibrations.  There  are  three  general  features  of  an 
alluvial  stream  that  a  calibrated  mathematical  model  should  reproduce  at  least 
approximately,  when  compared  with  actual  field  observations.  These  features 
are: 

(1)  Water-surface  profiles  for  a  range  of  flow  discharges, 

(2)  Sediment-transport  quantities,  and 

(•>)  General  bed-profile  trends. 

The  first  step  was  to  develop  a  rating  curve  that  would  give 
Manning's  n  of  the  main-channel  sections  in  the  study  reach  for  the  complete 
range  of  flow  discharges  expected  in  the  MR.  Using  the  daily-averaged 
quantities  for  the  DMR  flow  discharge,  plus  the  water-surface  elevation 
and  the  flow  discharge  at  Canton,  the  HEC-6  model  was  run  for  one-day 
time  intervals  using  a  wide  range  of  flow  discharges.  For  each  daily- 
averaged  discharge  used,  a  value  of  n  was  selected  so  that  the  computed 
upstream  water-surface  elevation  would  correctly  match  the  measured  value 
at  L&D  19.  A  linear  interpolation  of  the  results  using  the  least -squares 
method  was  made,  the  result  of  which  is  shown  in  figure  II-A-3.1,  and  the 
relation  is  expressed  by 

n  =  9.30  x  10-<”  Q  +  0.0101  lor  Q  118,000  cfs 

n  =  0.0213  for  Q  •  118,000  cfs 


(2.4) 

(2.5) 


ween  Manning  s  n  and  flow 
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Note  that  these  relations  were  used  to  describe  the  main  channel  only. 

A  constant  value  of  n  =  0.06  was  selected  for  use  in  the  model  for  all 
overbank  subsections. 

From  the  rating  curve  for  a  test  set  of  daily  discharges,  it 
was  found  that  the  maximum  error  introduced  in  the  upstream  water-surface 
elevation  at  Keokuk  was  of  the  order  ±0.5  ft.  Because  the  water-surface 
elevation  at  the  Keokuk  hydroelectric  plant  varies  ±0.5  ft  during  the 
day  due  to  varying  power  generating  demands,  the  rating  curve  was  judged 
to  be  adequate  and  was  employed  throughout  the  study. 

After  determination  of  the  prediction  for  n,  the  next  step  was 
to  reproduce  the  sediment -transport  quantities  obtained  from  field  measure¬ 
ments.  The  HEC-6  program  was  run  using  flow  conditions  typical  of  those 
found  during  the  1978  field  study. 

The  program  significantly  underestimated  the  sediment  discharges 
measured  in  the  field.  Consequently,  the  water-sediment-inflow  relation¬ 
ship  for  the  DMR,  which  was  questionable  from  the  start,  was  modified  as 
follows.  On  the  days  that  sediment-  and  water-discharge  measurements  were 
taken  on  the  DMR  during  the  1978  study,  another  cross  section  that  was 
measured  downstream  on  the  same  or  a  consecutive  day  was  chosen.  The  simu¬ 
lated  sediment-discharge  input  from  the  DMR  was  compared  with  the  sediment 
discharge  measured  in  the  field.  Modified  suspended-sediment  and  bed-load 
discharges  were  then  computed  by  multiplying  each  of  the  suspended-sediment 
and  bed-load  discharges  predicted  from  the  original  relationship  by  the  ratio 
of  the  field-measured  and  simulated  quantities.  This  calculation  yielded 
a  modified  set  of  water-  and  sediment -discharge  quantities,  through  which 
least-squares  fits  gave  the  following  altered  water-sediment-inflow 
formulas: 

Q  =  8.27  x  10-6  Q2'40  (2.6) 

s  -9  2  89 

=  1.87  x  10  Q  (2.7) 

The  program  was  run  again  utilizing  the  new  formulas,  and  the 
computed  total  sediment  loads,  Qp  =  +  Qg,  that  resulted  were  found  to 

be  in  much  better  agreement  with  the  field  measurements.  Table  II-A-3.1 
shows  the  original  and  modified  sets  of  sediment -in flow  relations  used 
for  the  DMR  in  'he  HEC-6  computation  of  sediment  discharge  at  specific 
locations  in  the  MR  study  reach .  The  difference  at  RM  346.97  was  expected 


because  the  model  cannot  adequately  simulate  the  significant  sediment 
contribution  from  a  side  channel  upstream  from  this  section.  In  other 
words,  the  model  cannot  account  for  effects  due  to  channel  bifurcation. 

In  addition,  the  model  underestimated  the  sediment  discharges  during  the 
two  low-flow  periods.  However,  the  MR  is  in  a  corparatively  static  or 
stable  condition  during  low  flows,  so  the  discrepancies  occur ing  during 
similar  periods  of  low  flow  are  not  considered  critical  to  the  model  per¬ 
formance. 


Table  II-A-3.1 

Comparison  of  Total  Sediment  Discharges  Using  Original 
and  Modified  Des  Moines  River  Rating  Curves 


River 

Mile 

Water 

Discharge 

Q(cfs) 

Field 

Measured 

Qp  (tons/day) 

Original 

Qp  (tons/day) 

Modified 

Qp  (tons/day) 

346.97 

159,183 

223,544 

126,997 

157,872 

350.33 

142,010 

200,534 

124,916 

210,273 

354.89 

107,869 

106,509 

61,125 

103,335 

351.87 

47,257 

5,533 

1,783 

1,816 

346.97 

42,727 

4,621 

563 

565 

Finally,  the  calibrated  model  should  reproduce  general  bed-profile 
trends.  A  4-month  run  from  May  through  August  1978  was  compiled  to  acquire 
a  rough  idea  of  how  well  the  model  results  would  approximate  the  1978  field 
study  results.  Figure  II-A-3.2  shows  the  simulated-bed  fluctuation  of  the 
monthly-averaged  flows  of  June  and  August  1978  compared  to  the  results  from 
the  two  different  sampling  periods  during  the  1978  field  study.  This 
run  employed  general  cross  sections  obtained  from  the  topographic  maps. 

The  detailed  cross  sections  from  the  1976  field  study  were  included  in 
the  model  during  later  runs.  The  model  stows  fairly  good  agreement  with 
field  measurements.  Hie  bed  elevation  from  field  measurements  was  approxi¬ 
mated  by  subtracting  the  mean  depth  (i.e.,  cross-sectional  area  divided  by 
river  width)  from  the  water-surface  elevation.  Hie  bed  elevation  of  the 
second  trip  (14-24  August  1978)  was  then  subtracted  from  the  bod  elevation 
of  the  first  trip  (26  June-11  July  1978)  at  the  resjxictive  study  sections 
to  obtain  the  points  plotted  in  the  figure.  Hie  figure  displays  the  fairly 


Figure  II-A-3.2  Comparison  of  computed  and  measured  bed-elevation  changes 
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erratic  changes  that  were  found  over  short  distances  during  the  field 
study.  The  differences  clearly  demonstrate  the  problem  that  is  encoun¬ 
tered  when  trying  to  compare  three-dimensional  phenomena  and  measurements 
acquired  from  field  studies  to  simulated  results  that  are  obtained  from 
one-dimensional  mathematical  models.  Note  that  after  the  model  calibration, 
additional  cross  sections  were  compiled  from  the  1976  field  study  to  in¬ 
clude  more  detail  with  respect  to  side  channels  and  cross-section  geometry 
in  the  model. 

In  suirmary,  the  calibration  procedure  entailed  first  formulating  a 
rating  curve  that  would  give  an  indication  of  the  main-channel  roughness 
through  a  value  of  Manning’s  n  for  a  given  flow  discharge.  The  resulting 
relation  produced  a  maximum  error  between  the  computed  and  daily-averaged 
upstream  boundary  water-surface  elevation  of  not  more  than  ±0.5  ft.  This 
error  was  felt  to  be  insignificant  since  varying  power  demands  throughout 
the  day  at  the  hydro-electric  plant  produced  fluctuations  of  ±0.5  ft  in 
the  water-surface  elevation. 

Next ,  sediment  discharges  computed  from  simulation  runs  were  found 
to  be  much  less  than  the  sediment  discharges  measured  in  the  field.  Con¬ 
sequently  the  water-sediment- inflow  formula  for  the  DMR,  which  was  initially 
questionable,  was  revised  to  better  approximate  the  field  measurements, 
while  tlx?  MR  water-sediment-inflow  relation  at  LfrD  19  was  left  unchanged 
since  it  was  based  on  a  10-vear  record  and  was  considered  to  be  accurate. 
Sediment  discharges  computed  from  simulation  runs  with  the  new  DMR  rela¬ 
tions  better  approximated  sediment  discharges  measured  in  the  MR,  although 
the  discharges  were  still  underestimated  during  low-flow  periods.  Because 
the  higher  flow  regimes  seemed  to  be  adequately  modeled  and  are  considered 
more  important  when  trying  to  represent  a  dynamic  river  system,  while  lower 
flow  regimes  were  relatively  stable  with  little  change  in  river  conditions, 
no  further  revisions  were  made. 

The  final  check  was  to  observe  the  general  trends  of  tlx?  bed-profile 
variation  with  time.  Results  from  a  simulation  run  .showed  scour  and  depo¬ 
sition  trends  similar  to  the  1978  field-study  observat ions .  The  final 
revisions  of  the?  calibv  xl  model  were  inclusion  of  additional  ;ind  more 
detailed  cross  sections  at  locations  identical  to  those  studied  during 
the  1976  study,  and  incorporation  of  a  dredging  operation  that  the  00K  (HI) 
iierformcxl  on  A  September  1978  near  tlx*  tip  of  Hazard  Island  Ixdwt'en  RM  !H9 
and  850. 
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It  should  be  mentioned  that  more  research  into  the  effects  of  varying 
the  bed-material  compositions,  the  percentages  of  the  various  size  frac¬ 
tions  composing  the  total  transported  sediment  load,  and  water-sediment - 
inflow  formulas  could  produce  a  better-calibrated  model. 

4.  Sensitivity  analysis. 

a.  Values  of  Manning's  n.  The  effect  of  values  of  Manning's  n  on 
the  water -surf ace  profile  can  be  seen  in  figure  II-A-4.1.  Varying  n  from 
0.0164  to  0.0210  changes  the  upstream  water-surface  elevation  by  only  1.4 
ft  for  a  flow  discharge  of  80,600  cfs  in  a  21-mi  reach.  However,  the  impor¬ 
tance  of  the  water-surface  elevation  is  in  its  involvement  in  the  calcula¬ 
tion  of  the  capacity  of  the  flow  to  transport  the  sand-sized  particles. 
Figure  II-A-4.2  shows,  for  identical  water  discharge,  how  much  the  sand- 
sized-sediment  discharge  may  vary  depending  on  the  value  of  n  used  in 
calculating  the  sediment-transport  rate.  The  silt-  and  clay-sized  sedi¬ 
ment  discharges  depend  almost  entirely  on  the  water  discharge,  and  so  are 
not  affected  by  n.  The  amount  of  sand  being  transported,  deposited, 
and  scoured  determines  the  consequent  thalweg-elevation  shift,  which  in 
turn  affects  the  sand-transport  capacity  of  the  flow.  Hence,  because 
of  the  intricate  dependence  among  the  variables  described  here,  the  rela¬ 
tionship  involving  n  should  be  formulated  from  reliable  field  data. 


b.  Water  temperatures .  Water  temperatures  are  used  by  the  computer 
program  in  the  calculation  of  fall  velocities.  The  emphasis  is  placed 
on  the  sand-sized  sediment  particles,  which  determine  the  resulting  bed 
profile.  Figure  II-A-4.3  exhibits  the  difference  in  the  transported  sand 
discharges  for  two  fairly  extreme  temperatures.  However,  the  water  tem¬ 
perature  seldom  undergoes  large  fluctuations  within  a  week,  so  rough  esti¬ 
mates  of  temperatures  may  be  used  with  reasonable  confidence. 


e.  P/Kl-raiteria]  profiles.  Several  short  simulation  runs  were 
made  using  different  bed-material  eornjx>si t ions  to  examine  the  effects  of 
file  lx ‘d-mafer i a  1  size  distribution  on  the  a xnput.at ional  results.  1’ 
signi  f  ieimLly  different  Ixd-maferial  profiles  were  pnxlueed  with  the 
median  size,  1)^,  ranging  fr<xn  0.35  urn  to  1.12  nm.  One  profile  was 
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Figure  II-A-4.1  Effect  of  Manning’s  n  on  water-surface  profiles 
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Figure  II-A-4,3  Effect  of  temperature  on  sand  load 
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inserted  into  five  cross  sections  covering  a  four-mile  reach  and  run  for 
a  seven-day  trial  period  using  a  flow  discharge  of  42,727  cfs.  The  other 
three  bed-material  profiles  were  then  alternately  inserted,  and  similar 
trial  runs  were  made.  The  results  obtained  shoved  that  the  coarser  bed- 
material  profile  (larger  D^q)  produced  greater  sand  discharges,  in  the 
vicinity  of  the  test  reach.  However,  the  differences  in  these  sand  dis¬ 
charges  were  not  significant,  and  the  thalweg  elevation  and  silt  and  clay 
discharges  remained  practically  constant  for  all  runs.  Note  that  a  low 
water  discharge  was  used  and  that  the  use  of  higher  water  discharges  would 
possibly  have  produced  mere  significant  differences.  In  the  example 
cited  above,  it  did  not  appear  that  the  program  was  particularly  sensi¬ 
tive  to  the  bed-material  compositions  initially  used  as  input  data.  This 
lack  of  sensitivity  is  advantageous  since  there  is  no  standard  procedure 
available  to  guide  selection  of  the  number  of  bed-material  samples  from  a 
field  study  necessary  to  compile  a  composite  size-distribution  profile 
for  the  entire  cross  section.  However,  the  ability  of  the  HEC-6  program 
to  incorporate  a  full  bed-material  profile  places  it  above  most  other  mathe¬ 
matical  models  that  generally  require  only  one  representative  size,  commonly 
taken  as 

The  HEC-6  program’s  ability  to  include  a  full  bed-material  pro¬ 
file  seems  to  permit  it  to  better  simulate  the  scouring  and  deposition 
process,  and  the  armoring  process.  Furthermore,  although  the  initial  bed- 
material  compositions  need  not  be  particularly  accurate,  the  final  profiles 
obtained  after  long-run  durations  would  indeed  embody  the  deposition  and 
scouring  effects  from  the  numerous  sediment -size  fractions  input  to  the 
model.  However,  as  will  be  seen,  the  characteristic  parameters  describing 
the  bed-material  profiles  can  vary  widely  over  short  periods  of  time 
and  over  short  distances.  Hence,  the  advantage  of  a  model’s  capability  to 
include  a  full  bed-material  profile  is  that  the  interaction  of  a  wide  range 
of  particle  sizes  will  be  represented  during  the  simulation  of  the  depo¬ 
sition  and  scouring  pnxx?sses.  Therefore?,  the  model  would  be  able  to  re¬ 
flect  actual  river  phenomena  to  a  more  realistic  extent. 

r>.  Discussion  of  results  from  the?  simulation  runs. 

a.  Bed  profiles.  Figure  ri-A-f>.l  si  lows  the  final  thalweg  ele¬ 
vations  obtained  from  the  runs  using  the  monthly-,  wtx'kly-,  ;md  daily- 
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Figure  II-A-5.1  Comparison  of  computed  thalweg  elevations  with  field 
data  (HEC-6  model) 
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averaged  flow  quantities  for  the  period  May  1976  through  August  1978.  The 
thalweg  is  defined  as  the  deepest  point  in  a  channel  cross  section.  In  gen¬ 
eral,  there  are  only  arall  differences  among  the  results  of  the  three  simu¬ 
lation  runs.  The  similarity  is  valuable  since  the  cost  of  running 
the  program  using  the  weekly-  and  daily-averaged  quantities  was  roughly 
two  and  eight  times  the  cost  of  the  monthly-averaged  run,  respectively. 

The  agreement  between  the  model  results  and  the  1978  field-study  results  is 
good  in  those  areas  where  the  rrcdel  was  constructed  from  the  cross  sections 
surveyed  during  the  1976  field  study,  but  poor  in  the  other  areas  where 
the  cross  sections  were  established  from  the  old  topographic  maps.  This 
outcome  stresses  the  importance  of  having  accurate  cross-section  geometries 
obtained  from  field  surveys,  from  which  the  initial  conditions  of  a  mathe¬ 
matical  model  can  be  determined. 

Figure  II-A-5.2  shows  the  net  change  between  the  final  and  initial 
bed  profiles  for  each  of  the  three  runs.  A  positive  change  denotes  depo¬ 
sition,  while  a  negative  change  denotes  scouring.  The  model's  overall  trends 
agree  well  with  the  observations  from  the  1976  and  1978  field  studies.  The 
areas  near  Fox  Island  (Rm  355)  and  Buzzard  Island  (EM  349.5)  are  known  for 
their  recurring  shoaling  problems,  which  are  indicated  by  the  model  re¬ 
sults,  although  the  7.6  ft  of  deposition  computed  using  the  monthly-averaged 
quantities  is  unreasonable.  Deposition  was  also  particularly  dominant 
at  RM  347  during  the  1978  field  study,  and  it  is  also  exhibited  by  the 
model  results.  Finally,  deposition  would  be  expected  to  occur  at  the  down¬ 
stream  boundary  of  the  model  at  L&D  20,  due  to  the  particle  settling  effect 
of  the  backwater  pool.  For  a  supplementary  comparison ,  the  differences 
in  the  thalweg  elevations  and  the  mean  bed  elevations  (i.e.,  water-surface 
elevation  minus  the  mean  flow  depth,  d  =  A/W)  between  the  1976  and  1978  field 
studies  are  shown  in  figure  II-A-5.3. 

b.  Water-surface  profiles.  Figure  II-A-5.4  shows  the  varia¬ 
tion  in  the  recorded  and  confuted  upstream  water-surface  elevations  at 
liSJ)  19  obtained  from  the  simulation  run  using  the  dni  ly-averagr*!  flow 
qiuint i t ies.  Note  that  only  one  |x>int  |x>r  vwn'k  has  Ixx'ii  plotted.  ’Hie  ngnx'- 
mont  is  gcxxi  except  near  days  580  to  700 ,  representing  I  Deccmlx’r  1977  to 
1  April  1978.  This  disagreement  is  believed  to  be  due  to  the  effect  of 
the  ice  cover  in  Pool  20,  which  prevented  the  water  surface  from  rising  to 
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its  free-surface  level.  Consequently,  a  larger  storage  volume  is  required 
in  the  upstream  reach  of  the  pool,  thus  causing  the  actua1  recorded  water 
surface  elevation  at  L&D  19  to  be  much  higher  than  would  be  normally  expect¬ 
ed  for  these  flow  conditions.  The  overall  agreement  demonstrates  the 
adequacy  of  the  n-Q  relation  that  was  formulated. 

Figure  II-A-5.5  presents  three  representative  water-surface  pro¬ 
files  for  a  high,  intermediate,  and  low  water  discharge  obtained  from  a 
simulation  run  using  daily-averaged  flow  quantities.  Similar  profiles 
wore  obtained  using  the  monthly-  and  weekly-averaged  flow  quantities.  As 
seen  in  the  figure,  the  form  of  a  dip  in  the  water-surface  profile  occured 
at  EM  356  for  all  water-surface  profiles.  The  exact  cause  of  this  dip  is 
not  known,  although  it  is  thought  to  be  caused  by  the  sudden  channel  con¬ 
traction  experienced  between  RM  357  and  356,  with  a  subsequent  channel  ex¬ 
pansion  tetween  RM  356  and  355.59.  The  dip  did  not  appear  prior  to  the 
inclusion  of  the  additional  cross  sections. 

c.  Trap  efficiencies,  computer  time,  and  dredging  volumes. 

Table  II-A-5.1  lists  for  each  of  the  three  runs  the  total  sediment  inflow 
and  outflow  volumes,  the  resulting  trap  efficiency,  and  the  computer  time 
required  for  the  simulation.  As  can  be  seen  by  comparing  the  three  runs, 
the  amount  of  inflowing  sediment  volumes  increased  when  shorter  time  steps 
were  used.  Although  the  sand  discharge  composed  only  ten  percent  of  the 
total  sediment  volume,  the  sand  fraction  had  an  influential  role  in 
determining  deposition  and  scour  of  the  river  bed.  Note  that  the  IffiC-6 
program  does  not  consider  the  removal  of  silt  and  clay  from  the  bed 
once  deposition  occurs. 

The  weekly-  and  daily-averaged  flow  runs  required  roughly  one- 
third  more  and  four  times,  respectively,  the  amount  of  computer  time  re¬ 
quired  by  the  monthly-averaged-flow  runs.  In  terms  of  final  cost,  however, 
the  weekly  and  daily  runs  were  approximately  two  and  eight  times  the  cost 
of  the  monthly  runs,  respectively.  However,  the  results  from  the  three 
different  runs  do  not  vary  significantly.  Therefore,  to  save  time  and 
money,  the  monthly-  and  weekly-averaged  flow  quantities  should  suffice 
in  simulating  scour  and  deposition  trends,  al  tfmugh  during  jioak  flows  if 
is  advisable  to  include  several  periods  of  dai ly-averagod  flow  quantities. 


365  360  355  350  345  343 
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re  Il-A-5.5  Representative  water -surf ace  profiles  obtained  using  daily-averaged  flow  quantities 
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d.  Bed-material  profiles.  The  bed-material  profiles  are 
probably  the  most  difficult  river  characteristic  to  simulate  and  to 
compare  with  profiles  derived  from  field  studies.  The  size  distribu¬ 
tion  of  bed  material  can  vary  significantly  across  a  section,  where  depo¬ 
sition  may  be  occurring  in  one  area,  while  another  area  may  be  undergoing 
scour.  Numerous  samples  were  taken  at  a  cross  section  in  an  attempt  to 
obtain  a  representative  composite  bed-material  profile.  The  HEC-6  program 
is  able  to  incorporate  a  full  bed-material  profile  at  each  cross  section. 
However,  table  II-A-5.2  shows  how  much  variation  can  be  expected  in  the 
fundamental  bed-material  parameters  over  time  and  distance.  The  data 
representing  the  three  sampling  trips  made  during  the  1970  field  study 
indicate  the  variations  that  can  be  expected  due  to:  (1)  the  changes  that 
the  bed-material  composition  undergoes  with  time  (in  this  case,  the  four- 
month  period  from  May  to  September ,  1976)  and,  (2)  the  changes  resulting 
from  the  different  locations  of  sampling  sites  at  a  particular  channel 
cross  section.  The  columns  labeled  HEC-6  denote  the  resulting  parameters 
of  the  three  28-month  simulation  runs.  The  last  two  columns  represent 
the  results  of  the  1978  field  study.  Ideally,  one  would  expect  the  HEC-6 
results  to  be  similar  to  the  1978  field-study  results,  although  the 
slightly  different  cross-section  locations  could  cause  a  difference. 

Overall,  these  results  demonstrate  the  difficulty  in  modeling  and  comparing 
the  various  bed-material  parameters,  which  have  been  shown  to  vary  with  time 
and  location. 

One  final  point  worth  mentioning  refers  to  differences  obtained 
from  the  use  of  different  computer  systems  and  languages.  Running  an 
identical  sample  program  supplied  by  the  HEC  on  an  IBM  370  system  yielded 
slightly  different  results  than  those  obtained  by  the  HEC.  Also,  the  use 
of  standard  FORTRAN  and  Minnesota  FORTRAN  gave  slightly  different  results. 
However,  the  variations  found  were  of  the  order  of  t  2  percent,  and  thus 
are  not  considered  to  be  significant. 

B.  The  CSU  Mathematical  Models. 

1.  General  model  description.  Two  mathematical  simulation 
models,  developed  at  Colorado  State  University  (CSU),  Fort  Collins,  Color¬ 
ado  ((hen  and  Simons,  1980),  have  lx*en  appl  ied  to  the  l>x)1  20  st  udy  reach 
a.  I,  the  request  of  the  l  1 1  IK  l.o  invest  igale  I  h«  *  elTect  i  veness  el  I  he  models 


L  k 


349.24 

348.98 


99 


in  simulating  flow  conditions  and  phenomena.  The  first  model,  a  one¬ 
dimensional  uncoupled-unsteady  water  and  sediment  routing  model  (UUWSR 
Model)  employs  an  implicit  numerical  method  for  water  routing  to  solve 
the  water  continuity  and  momentum  equations  assuming  a  fixed  bed.  The 
sediment  continuity  equation  for  sediment  routing  is  then  solved  at  the 
same  time  step.  The  UUWSR  model  lias  been  applied  previously  to  Pools 
4  through  8  in  the  Upper  MR  system  (Simons  et  al. ,  1979).  The  second 
model,  a  one-dimensional  steady-uncoupled  sediment  routing  model  (SUSR 
Model)  also  assumes  a  fixed  bed,  then  computes  the  backwater  profile  for 
a  step  discharge  by  solving  the  energy  equation.  The  bed-elevation  changes 
are  determined  at  the  end  of  the  time  step  by  solving  the  sediment  con¬ 
tinuity  equation.  The  SUSR  model  has  been  applied  to  the  Yazoo  River 
Basin  (Simons  et  al.  1978)  and  has  been  found  to  be  excellent  in  studying 
long-term  changes  in  a  complex  river  system. 

As  seen  with  the  HEC-6  model,  the  first  objective  was  to  use 
historical  records  in  addition  to  data  collected  from  the  1976  field  study 
(Nakato  and  Kennedy,  1977)  to  establish  the  initial  conditions  for  the 
models.  Simulation  runs  were  then  carried  out  through  the  time  of  comple¬ 
tion  of  the  1978  field  study  (Vadnal,  1979),  thereby  allowing  a  comparison 
between  the  field  results  and  those  obtained  from  the  model  predictions. 

Section  2  reports  on  the  construction  of  the  two  models  using 
the  available  geometric,  sediment,  and  hydrologic  data,  while  section  3 
discusses  the  calibration  of  the  models.  Section  4  analyzes  the  sensitivi¬ 
ties  of  the  two  models,  and  section  5  evaluates  the  performance  of  the  two 
models  and  the  results  obtained. 

2.  Construction  of  the  models. 

a.  Data  availability.  Field  data  from  the  two  field  studies 
conducted  by  the  IIHR  wore  used  along  with  historical  records  to  compile 
the  input  conditions  for  the  models.  Geometric  data  were  obtained  from 
the  1977  report  for  cross  sections  in  the  problem  shoaling  area,  while 
the  other  model  cross  sections  wore  obtained  from  a  1945  survey  conducted 
by  OOK(RI).  Water  ;uid  sediment,  inflows  Cor  the  DMR  wore  ohlaintxl  from 
records  ;uul  mexsuremenls  at  .St.  I'ranc.isvi lie,  Missouri.  Hydrologic  data 
(water  discharge  and  stage  elevations)  for  IjSJI  ID  and  L&D  20  wore  com¬ 
piled  from  discharge  and  stage  hydrographs .  Sediment-discharge  relationships 
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at  L&D  19  were  obtained  from  measurements  taken  there.  A  preliminary 
relationship  was  derived  for  the  study  reach  using  the  previously  discussed 
HEC-6  program  (U.S.  Army  Corps  of  Engineers,  1977).  Details  of  dredging 
operations  performed  during  the  model  study  period  (May  1976  through  August 
1978)  also  were  available. 


b.  Model  construction.  The  Pool  20  river  reach  was  divided 
into  27  cross  sections  with  space  intervals  ranging  from  0.2  mi  in  the 
Fox  Island  shoaling  area  to  1.2  mi  at  L&D  19.  Figure  II-A-1.1  displays 
the  locations  of  the  chosen  cross  sections  and  a  map  of  the  study  area. 

Since  the  wash  loads  typical  in  the  MR  are  generally  considered 
to  be  transported  continually  throughout  the  main  flow,  only  the  bed- 
material  discharges  affect  the  river-bed-elevation  changes.  Consequently, 
only  the  bed-material  discharges  were  considered  in  the  models  for  sediment 
routing.  However,  the  1976  field  study  found  that  roughly  10  percent  of  the 
suspended-load  samples  collected  near  the  shoaling  areas  was  in  the  sand-size 
range,  while  roughly  1  percent  of  the  suspended- load  samples  collected 
near  L&D  19  was  in  the  sand-size  range,  and  approximately  25  percent  of 
similar  samples  taken  in  the  DMR  contained  sand-size  particles.  Thus, 
bed-material-transport  relationships  were  derived  from  the  available  data 
which  incorporated  the  appropriate  fraction  of  the  sand-sized  particles 
found  in  suspension.  The  relation  for  the  general  study  reach  is 

gmI  =  9.82  x  10"6  V3-58  (2.8) 


where  g^  is  the  unit  width  bed-material  discharge  in  cfs/ft,  and  V  is 
the  mean  flow  velocity  in  ft/s.  This  equation  was  used  in  both  the  UUWSR 
and  SUSR  models  to  obtain  the  sediment-transport  rate  in  Pool  20  by  multi¬ 
plying  ggjj  by  the  top  channel  width. 

For  the  sediment  input  from  L&D  19, 


%I(LD19) 


=  8.42  x  10  7  V3’5 


(2.9) 


\vh* nv  is  the  unit  width  Ixxl-mcitorial  discharge  at  L&D  19  in 

MM(  I J  )1 9) 

cts/fl  ,  and  V  is  tlu*  mean  I'low  Y<  lo<  ily  in  ll/s. 
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For  the  sediment  input  from  the 


=  5.7  x  10"8  Q1,74 


(2.10) 


where  is  the  bed-material  discharge  in  cfs,  and  Q  is  the  flow 

discharge  in  cfs  for  the  DMR. 

The  preliminary  n-Q  relationship  for  Pool  20  obtained  by  the 
HEC-6  pregram  was  given  previously  as: 


n  =  9.502  x  10  Q  +  0.01  for  Q  <  118,000  cfs 


n  =  0.0213 


for  Q  >  118,000  cfs 


(2.4) 


(2.5) 


However,  the  UUWSR  and  SUSR  models  require  a  power  relation  of  the  form: 


n  =  nQ(aQ  ) 


(2.11) 


Thus,  (2.4)  was  approximated  by 


n  =  0.0213  (0.0324)Q 


(2.12) 


At  the  beginning  of  the  calibration  procedure  (2.5)  and  (2.12)  were  used, 
but  were  later  modified  to  obtain  better  agreement  between  measured  and 
computed  stage  elevations  at  L&D  19. 

3.  Calibration  of  the  models.  A  calibrated  mathematical  model 
should  be  able  to  simulate  flow  characteristics  and  geomorphic  changes  in 
the  modeled  study  reach.  Three  important  factors  for  simulation  are: 


(1)  The  water  discharge  and  water-surface  elevation  at  computa¬ 
tional  cross  sections, 

(2)  The  cross-sectional  changes,  and 

(3)  The  sediment-transport  rates. 
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a.  Calibration  of  the  UUWSR  model.  The  relations  involving  n 
were  modified  to  obtain  better  agreement  between  the  measured  and  computed 
stage  hydrographs  at  the  upstream  boundary,  L&D  19,  from  1  May  1976  to 

31  August  1976.  The  modifications  yielded: 

n  =  0.0213  (0.211)Q°‘ 13  for  Q  <  180,000  cfs  (2.13) 

n  =  0.0219  for  Q  >  180,000  cfs  (2.14) 

The  utilization  of  these  relations  gave  differences  in  stage  elevation 
at  the  upstream  boundary  typically  less  than  1  ft,  except  during  the  period 
from  December  1977  to  March  1978  where  ice-jam  effects  caused  larger  dif¬ 
ferences. 

Cross-sectional  changes  were  computed  at  one  Fox  Island  location 
and  at  four  locations  near  Buzzard  Island.  The  model  was  calibrated  to 
reproduce  the  geomorphic  changes  at  these  locations.  The  results  of  the 
computations  at  the  five  locations  are  shown  in  figures  II-B-3.1  through 
II-B-3.5  together  with  the  cross-sectional  shapes  found  at  the  completion 
of  tiie  field  study  in  August  1978.  Note  that  for  this  one-dimensional  model, 
the  overall  bed-area  changes  were  computed  at  each  cross  section,  and  then 
distributed  according  to  relative  conveyance  to  compute  the  new  cross-section 
shapes  shown  in  the  figures.  This  sediment -distribution  method  is  de¬ 
scribed  by  Simons  and  Chen  (1979). 

Modifications  of  the  sediment-transport  relations,  (2.3);  (2.9); 
and  (2.10),  to  produce  a  better  calibrated  modtl  were  not  attempted  due  to 
the  lack  of  more  recent  detailed  information  regarding  cross-section 
shapes,  and  the  fact  that  the  1976  and  1973  field  studies  surveyed  dif¬ 
ferent  cross  sections.  However,  correct  trends  in  the  thalweg  elevations 
were  predicted  by  the  model  in  the  recurrent  shoaling  areas  near  Fox  Island 
and  Buzzard  Island,  as  can  be  seen  in  figure  II-B-3.6.  Thus,  it  appears 
that  the  calibrated  model  drx)«  adequately  simulate  geomorphic  changes  over 
the  study  reach. 

b.  Calibration  <>l  tlx*  SUSP  nixlel  .  A  similar  pnxxKlure  In  that 
(liwriU'd  lor  I  he  IH1WS1J  Model  was  used  in  I  he  ni  >d  i  !  im  I  ions  ol  I  lx  *  relations 


ISTANCE  FROM  LEFT  BANK  (FT) 

omputed  cross-section  profiles  at  RM  354.93 


DISTANCE  FROM  LEFT  BANK  (FT) 

computed  cross-section  profiles  at  RM  349.82 


DISTANCE  FROM  LEFT  BANK  (FT) 

Measured  and  computed  cross-section  profiles  at  RM  349.45 


STANCE  FROM  LEFT  BANK  CFT) 

nnutcd  cross-sect  ion  profiles  at  RM  349.29 


STANCE  FROM  LEFT  BANK  (FT) 

nputed  cross-section  profiles  at  RM  348.96 
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involving  n.  The  revised  formulas  are  given  by 

n  =  0.0213  (0.211)  Q0-14  for  Q£  133,000  cfs  (2.15) 

n  =  0.0234  for  Q  _>  133,000  cfs  (2.16) 

These  values  calculated  for  n  using  (2.15)  and  (2.16)  are 
about  7  percent  larger  than  those  values  calculated  using  (2.13)  and  (2.14) 
in  the  UUWSR  model.  The  difference  is  due  to: 

(1)  The  SUSR  model  uses  an  energy  equation  for  backwater  compu¬ 
tations,  while  the  unsteady  waterfront inuity  and  momentum  equations  are 
used  in  the  UUWSR  model. 

(2)  The  SUSR  model  approximates  the  geometric  properties  at  a 
computational  cross  section  using  a  power  relation,  while  the  UUWSR  model 
uses  linear  interpolation  of  stages  versus  geometric  variables  (cros^- 
sectional  area,  top  width,  and  conveyance). 

As  with  the  UUWSR  model,  differences  between  the  recorded  and 
computed  stage  elevations  at  L&D  19  for  the  study  period  were  typically 
less  than  1  ft,  although  ice- jam  effects  caused  larger  differences  between 
December  1977  and  March  197S. 

Cross-sectional  shapes  that  were  computed  using  the  SUSR  model 
are  also  displayed  in  figures  II-B-3.1  through  II-B-3.5,  and  they  were  com¬ 
pared  to  the  final  shapes  determined  during  the  1973  field  study.  As 
discussed  earlier,  no  further  alterations  were  made  to  the  sediment -transport 
equations  to  obtain  a  better  calibrated  model.  However,  figure  II-B-3.7 
shows  that  the  SUSR  model  also  correctly  predicted  general  aggradation 
trends  in  the  vicinity  of  the  recurrent  shoaling  area  near  Fox  Island. 

Hence,  it  appears  that  the  SUSR  model  also  may  be  used  to  simulate  geomor- 
phic  changes  over  the  study  reach. 

4.  Sensitivity  analysis. 

a.  Spatial  design.  Spatial  desij^n  refers  to  a  model's  repre¬ 
sentation  ol’  the  physical  characteristics  of’  a  river  system,  including 
such  in  format  ion  as  the  location  of  tributaries  and  data  on  channel 
properties.  In  sediment  studies,  the  spacing  of  cross  sections  can  bo 

estimated  in  terms  of  a  particle  settling  lon'd-h  which  is  given  by 
\) 

w 


v 


(3.17) 


eld  data  (SUSR  model) 
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where  x  is  the  particle  settling  length;  D  is  the  depth  of  flow;  V  is 
the  velocity  of  flow;  and  w  is  the  fall  velocity  of  the  median  bed-material 
sediment  size,  D^q.  Since  the  largest  amounts  of  sediment  are  transported 
during  high  flows,  typical  high-flow  values  should  be  used  in  (2.17).  Using 
a  recaimended  spacing  from  lx  to  lOx  (Simons  et  al . ,  1979a),  the  spacing 
of  0.2  mi  to  1.2  mi  employed  in  the  models  was  found  to  be  adequate,  although 
effects  from  varying  the  cross-section  spacing  were  not  investigated.  Note, 
however,  that  cross  sections  should  be  located  in  areas  of  interest,  at 
control  points,  and  at  places  of  sudden  changes  in  water-surface  profiles. 

b.  Temporal  design  and  operational  discharge.  Temporal  design 
refers  to  a  model's  ability  to  simulate  changes  in  the  water  and  sediment 
inputs  to  a  river  system  with  time,  as  well  as  representing  the  changes  in 
the  water  and  sediment  discharges  that  occur  throughout  the  river  system. 
Temporal  variations  are  generally  mild  in  the  Pool  20  study  reach  of  the 
MR,  and  since  discharge  hydrographs  are  given  on  a  daily  basis,  the  use  of 
a  basic  one-day  time  step  was  considered  sufficient  and  convenient  for 
the  model.  However,  because  changes  occurring  along  the  river  bed  are  mini¬ 
mal  during  low  flows  when  sediment-transport  rates  are  small,  larger  time 
steps  may  be  used  or  smaller  discharges  can  be  bypassed  to  increase  the 
effectiveness  of  the  model.  As  a  result  of  an  analysis  performed  on  the 
absolute  values  of  the  changes  in  thalweg  elevations  with  various  discharges , 
those  portions  of  the  discharge  hydrograph  having  flows  less  than  40,000  cfs 
were  bypassed. 


c.  Effects  of  time  steps.  The  average  of  the  absolute  values 
of  the  differences  between  the  computed  and  recorded  stage  elevations  at 
L&D  19  were  compared  for  four  different  time  intervals  (1  day,  3  days,  5 
days,  and  10  days).  The  results  indicated  that  the  maximum  time  step  that 
can  be  used  in  the  UUWSR  model  witliout  losing  significant  accuracy  is  about 
5  days.  Similar  results  were  obtained  when  the  absolute  cumulative  changes 
of  the  computed  thalweg  elevations  at  each  cross  section  versus  time  were 
compared  for  the  different  time  intervals. 

Similar  comparisons  wore  made  using  the  SUSR  model,  liut  with 
time  intervals  of  1  day,  r>  days,  10  fin.ys,  and  30  days.  (VmpuU'd  results 


were  poor  when  the  30-day  time  interval  was  used.  However,  the  utilization 
of  a  variable  time  step  enploying  longer  intervals  for  smaller  discharges  and 
shorter  intervals  for  larger  discharges  produced  good  agreement  with  the 
results  from  a  1-day  time  step.  Hence,  it  was  shown  that  the  SUSR  model 
could  adequately  use  a  variable  time  step  and  still  produce  accurate  results. 
It  should  be  noted  that  the  SUSR  model  required  less  computer  time  than  the 
UUWSR  model. 


5.  Discussion  of  results.  Doth  the  UUWSR  and  SUSR  models  show 
promise  in  their  ability  to  predict  river  changes  in  Pool  20  of  the  Upper 
MR.  The  UUWSR  model  provided  a  better  stage-elevation  prediction,  and  both 
models  calculated  comparable  bed-elevation  changes.  The  SUSR  model  was  able 
to  use  a  longer  time  step  and  required  less  computer  time,  implying  that  the 
SUSR  model  is  better  suited  than  the  UUWSR  model  for  studying  long-term 
inpacts. 

C.  The  CHAR2  Mathematical  Model. 


1,  General  model  description.  A  one-dimensional  mathematical 
model,  CHAR2,  was  developed  by  Sogreah,  a  consulting  firm  in  Grenoble, 
Trance  (Cunge  and  Perdreau,  1972;  Sogreah,  1981).  The  model  consists  of 
a  one-dimensional  steady-flow  equation  and  a  sediment -continuity  equation: 


(%  +  gy)  +  =  o 

3x  dA2  D r 

(1  ■  p)  H +  h  !£ =  0>  G  =  G(y>z’b>d- 
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(2.18) 

(2.19) 


where  Q  =  Q(x)  =  flow  discharge;  A  =  A(x)  =  cross-sectional  area;  y  -  y(x,t) 

=  water-surface  elevation;  z  =  z(x,t)  =  river-bed  elevation;  D  -  conveyance 
factor;  p  =  porosity  of  river-bed  sediment;  b  =  river-bed  width  associated 
with  a  movable  bed;  G  =  G(x,t)  =  sediment  discharge;  d  =  sediment  character- 
istics;  and  g  -  gravitational  acceleration.  It  is  assumed  that  flow  celer¬ 
ities  are  much  greater  than  bed- form  movement,  and  the  flow  resistance 
is  given  by  the  Manning  roughness  coefficient  which  is  considered  constant 
in  time'.  The  bed  material  is  assumed  to  be  homogenous.  Ttie  model  considers 
only  Ixxl-load  tr;uisfx>rt,  allowing  the  use  of  one  of  t .ho  Moyor-IV'ter,  Ilngelund- 
tlnnseu,  Dubovs,  or  Kinstoni-Drown  formulas.  rlbe  numerical  scheme  utilized 


to  solve  the  system  of  nonlinear  hyperbolic  partial  differential  equations 
is  an  inf.  • icit  finite-difference  method  with  a  double-sweep  procedure 
developed  by  Cunge  and  Perdreau  (1972). 

2.  CHAR2  model  construction.  Application  of  the  CHAR2  nrxiel 
for  Pool  20  of  the  MR  first  required  a  general  model  layout  which  is  shown 
in  figure  II-C-2.1.  Hie  reach  between  RM  348  and  RM  357,  which  includes 
two  major  shoaling  area,  used  a  spatial  interval  of  0.5  mile,  while  the 
spatial  interval  for  the  rest  of  the  reach  ranged  from  0.8  mile  to  1.5 
mile.  The  input  data  and  boundary  conditions  utilized  were  as  follows: 

(1)  Input  data.  At  each  computational  point,  the  main-channel 
cross-section  geometry  was  constructed  such  that  the  cross-section  profile 
varied  laterally  in  a  stepwise  manner.  Note  that  only  main-channel  cross 
sections  were  used  in  the  model  construction.  At  each  subsection,  median 
bed-material  size  and  bed  elevation  were  given  using  available  data.  The 

porosity  of  river-bed  sediment  (0.49)  and  the  kinematic  viscosity  of  water 

-5  2 

(1.08  x  10  ft  /s)  were  assumed  to  be  constant  throughout  the  model  cali¬ 
bration  period.  The  constant  value  of  n  =  0.023  ft  -1^3s  was  adopted;  this 
value  was  based  on  the  calibration  result  of  the  HEC-6  program.  The 
Einstein-Brown  formula  was  used  to  model  bed  load. 

(2)  Boundary  conditions.  The  following  boundary  conditions  were 
given  for  the  model  construction: 

a.  At  L&D  19  (upstream  boundary),  the  monthly-averaged  MR  dis¬ 
charge,  Q,  was  given  as  a  function  of  time,  t,  for  the  model  calibration 
period.  Since  L&D  19  trapped  most  of  the  coarse  sediments,  no  sediment 
input  at  this  upstream  boundary  was  considered  (note  that  the  CHAR2  program 
considers  only  bed-load  transport). 

b.  At  the  MIR  mouth,  the  monthly-averaged  DMR  discharge  was 
given.  Total  sediment  input,  Qi ,  from  the  DMR  was  calculated  as  a  sum  of 
partial  suspended-load  discharge,  Q  ,  and  bed-load  discliarge,  Q^: 


<ir  “  %  0.25  Qs  (2.30) 

in  which  the  factor  of  0.25  was  employed  since  approximately  25  percent  of 
the  suspended  load  in  tin'  DMR  contains  sand  materials.  Hie  qiuint.i ties  C) 
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and  Qg  were  determined  as  functions  of  0  based  on  the  USGS  and  the  HUE 
field  data: 

Qg  =  3.10  x  10-3  Q1'66  (2.21) 

and 

=  1.19  x  10-8  Q2,59  (2.22) 

where  Q  is  in  units  of  cfs,  and  Qg  and  are  in  units  of  tons/day. 

c.  At  L&D  20  (downstream  boundary),  the  monthly-averaged  pool 
elevation,  y,  was  given  as  a  function  of  t. 

3.  Results  of  the  CHAR2  model  calibration.  The  computer  simula¬ 
tion  was  conducted  by  Sogreah  for  the  time  period  from  1  May  1976  to  31 
August  1978  (28  months).  As  described  in  the  preceding  section,  monthly- 
averaged  input  data  were  employed.  The  computer  run  utilized  a  temporal 
computation  interval  ranging  between  6  hrs.  and  5  days.  The  computer  print¬ 
out  provided  information  on  water-surface  and  thalweg  elevations,  mean 
velocity,  flow  discharge,  and  sediment-transport  rate  at  each  cross  section 
of  the  study  reach.  Figure  II-C-3.1  illustrates  the  initial  longitudinal 
thalweg  profile  and  calculated  water-surface  profile.  The  observed  water- 
surface  elevation  at  the  upstream  boundary  (RM  364.2)  agrees  very  well  with 
the  calculated  value.  Note  that  the  water-surface  elevation  at  the  downstream 
end  was  given  as  a  boundary  condition.  Figures  II-C-3.2  and  II-C-3.3  show 
similar  plots  for  t  =  185  days  and  t  =  543  days,  respectively.  Since  there 
are  no  field  data  available  on  thalweg  elevations  for  these  days,  it  is  diffi¬ 
cult  to  judge  the  accuracy  of  the  longitudinal  bed-profile  prediction.  How¬ 
ever,  the  overall  prediction  of  the  water-surface  profiles  seems  to  be  quite 
satisfactory.  Figure  II-C-3.4  provides  both  the  initial  and  calculated 
(t  =  850  days)  thalweg  elevations  at  all  computation  points  as  well  as 
those  which  were  measured  by  IIIIR  and  00E.  Although  the  model  considered 
only  the  bed- load- transport  rate,  the  overall  prediction  of  the  thalweg 
elevation  seems  to  be  satisfactory  except  at  certain  computational  points. 

A  rough  comparison  was  made  between  the  MR  lxxl-1oad  diseliar^s  computed  from 
the  (1IAR2  prognmi  and  Lite  approximately  cx|ui valenl.  lied- load  discharges  from 
the  l‘)7S  field  .study,  computed  l>y  si  mini  nr;  the  incnsunxf  lx  *d- load  discharg*  « 
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and  25  percent  of  the  measured  suspended-load  discharges.  The  comparison 
is  shown  in  table  II-C-3.1.  The  1978  study  cross  sections  chosen  were 
those  closest  to  the  cross  sections  modeled  by  CHAR2.  The  CHAR2  values, 
although  reasonable,  generally  underestimated  the  field  study  values.  How¬ 
ever,  the  accuracy  of  the  comparison  when  adding  25  percent  of  the  suspended- 
load  discharge  to  the  bed-load  discharge  is  unknown. 


Table  II-C-3.1 


Comparison  of  Bed-Load  Discharges  between  the  CHAR2 
Program  and  1978  Field  Study  at  Selected  Locations 


River 

Mile 

(RM) 

Date 

Q 

(cfs) 

V.25Qs 

(tons/day) 

QgCCHAI^) 

(tons/day) 

Q(CHAR2) 

=  148,842  cfs  Date 

:  29  June  1978 

355.0 

26  June 

107,869 

24,524 

20,148 

350.0 

30  June 

135,596 

40,024 

34,573 

349.5 

5  July 

101,830 

24,325 

29,765 

349.0 

7  July 

96,800 

16,361 

27,704 

347.0 

11  July 

159,183 

51,077 

28,391 

Q(QIAR2)  = 

=  62,004  cfs  Date: 

29  August  1978 

355.0 

16  Aug 

48,687 

1,115 

1,603 

350.0 

17  Aug 

38,077 

696 

458 

349.5 

21  Aug 

44,509 

3,394 

687 

349.0 

22  Aug 

30,703 

1,072 

687 

347.0 

24  Aug 

42,727 

1,050 

458 

D.  Sumnary  and  Recorrmendat ions .  The  present  study  consisted  of 
two  principal  parts.  PHASE  A  of  the  investigation  included  two  campaigns 
to  collect  and  analyze  detailed  field  data  in  the  Pool  20  study  reach  of 
the  Mississippi  River  (MR)  between  RM  347  and  355  and  in  the  Des  Moines 
River  (DMR)  near  its  mouth.  These  data  were  used  to  make  diagnoses  of 
the  recurrent  shoaling  problems  in  the  vicinities  of  Fox  Island  and  Buzzard 
Island.  The  empirical  flow- sediment -discharge  relationships  formulated 
using  the  data  enabled  the  determination  of  approximate  increases  in  sediment- 
transport  capacity  when  side  channels  were  closed  to  abate  shoaling  activities 
of  the  main  channel.  The  data  were  also  utilized  to  calibrate  various  numeri¬ 
cal  models  in  PHASE  B  of  the  study. 

PHASE  B  included  testing  of  the  HEC-G  model  (Hydrologic  Engineering 
Center,  Corps  of  Engineers),  the  UUWSR  and  SUSR  models  (Colorado  State 
University),  and  the  CHAR2  model  (Sogreah) .  The  HEC-6  program  was  run  at 
The  University  of  Iowa,  and  the  other  three  models  were  run  by  the  developers 
using  the  basic  initial  and  boundary  input  data  model  constructed  by  the 
Iowa  Institute  of  Hydraulic  Research  (IIHR).  Although  each  one-dimensional 
model  has  its  own  numerical  model  characteristics,  accurate  prediction  of 
a  longitudinal  river-bed  profile  required  them  to  have  in  common  the  fol¬ 
lowing  four  major  factors:  (1)  accurate  initial  conditions,  including  a 
cross-section  profile  and  bed-material  size  distributions  at  each  computa¬ 
tional  cross  section;  (2)  accurate  boundary  conditions  such  as  water  and 
sediment  inflows  along  the  model  boundaries,  quantitative  expressions  of 
suspended-  and  bed-load  discharges  and  information  on  the  size  of  the  sedi¬ 
ment  inputs,  and  stage  hydrographs  at  the  upstream  and  downstream  boundaries: 
(3)  bed-roughness  characteristics  at  each  computational  point;  and  (4) 
reliable  sediment-transport  formulas  which  describe  the  sediment-transport 
characteristics  in  the  study  reach.  It  is  extremely  important  to  understand 
the  interrelationship  between  these  factors;  an  accurate  estimate  of  sediment- 
transport  rate  depends  entirely  on  accurate  estimates  of  river-flow  charac¬ 
teristics  which  rcHiuuv  detain'd  geometric  information  as  well  as  nn  inlor- 
mat  ion  concerning  U'd  roughness,  which  adjusts  itself  according  to  the 
sedimont-tnuisport  rate.  'Hie  interaction  betwixn  the  flow  and  the  movable 
river  bed  is  a  continual,  dynamic  activity.  Therefore,  the  exclusion  of 
oven  one  item  listed  above  can  load  to  serious  errors  in  computer  simula¬ 
tions.  However,  since  one  can  hardly  be  provided  with  a  complete  set  of 
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input  data  in  a  practical  numerical  application,  a  number  of  assumptions 
often  have  to  be  made  to  close  the  gap  in  the  input  information. 

Unfortunately,  the  study  reach  lacked  certain  input  data  in 
varying  degrees;  the  most  serious  one  was  a  lack  of  information  on  geometric 
configurations  of  the  initial  MR  bed  profiles,  sediment-inflow  rates  from 
the  DMR,  and  bed-material  size  distributions  along  the  river. 

Simulation  runs  of  the  aforementioned  models  all  were  made  for  a 
28-month  time  period  between  May  1976  and  August  1978.  The  initial 
longitudinal  river-bed  profiles  for  the  HEC-6,  UUSWR,  and  SUSR  models  were 
constructed  mainly  using  CDE’s  1945  topographic  maps,  except  for  the  cross 
sections  measured  in  1976  by  I IIEFL  In  the  CHAR2  model,  more  recent  topographic 
data,  obtained  by  ODE  in  1974  and  1976,  were  incorporated  for  several 
sections.  Therefore,  the  initial  conditions  for  the  CHAR2  model  and  the 
other  models  were  slightly  different.  The  predicted  thalweg  elevations 
by  the  four  models  were  compared  with  the  measured  1978  values.  The 
degree  of  agreement  between  the  conputed  and  measured  values  seems  to  be 
almost  of  the  same  order  for  each  model.  Better  agreement  was  generally 
found  in  the  areas  with  sufficient  input  data.  As  far  as  the  averaging 
period  for  the  input  data  was  concerned,  monthly-averaged  input  data  were 
sufficient  in  both  the  HEC-6  and  CHAR2  models;  whereas  the  two  CSU  models 
required  a  5-day  time  step  for  a  flow  discharge  over  100,000  cfs,  a  10-day 
time  step  for  a  discharge  between  50,000  cfs  and  100,000  cfs,  and  a  30-day 
time  step  for  a  discharge  below  50,000  cfs. 

The  application  of  these  numerical  models  to  predict  and  evaluate 
more  accurately  the  river-bed  changes  in  the  study  reach  requires  the 
establishment  of  initial  bed  profiles  at  all  computation  points.  This  task 
can  be  accomplished  easily  by  detailed  soundings  including  side  channels 
along  the  reach.  Concerning  the  sediment -input  information  in  the  DMR, 
sediment  sampling  should  be  continued  at  St.  Francisville,  Missouri,  to 
establish  a  meaningful  and  reliable  flow-sediment  rating  curve  since  the 
UMR  is  the  major  source  of  sediment  responsible  for  the  recurrent  shoal  ing. 

With  these  sinple,  suppl  orient  ary  data  the  calibration  of  otic-dimensiona  1 
mrxiels  will  certainly  iH^oome  more  reliable,  and  I  he  long-f  enn  e!  I'eel  <»i 
side  channel  closures  roconinondt  ?d  i  n  Chapter  f  can  Ik*  losCd.  AMUni<h  i 
twr>~dimensional  mrxlel  has  Vx'on  recently  developed  by  (MU  and  tested  for 
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Fbol  4  of  the  MR,  the  future  of  such  models  is  still  in  the  dark  because 
of  the  lack  of  sufficient  field  data  to  calibrate  (note  here  that  there  are 
not  sufficient  input  data  for  even  a  QNE-D I  MENS  IONAL  model )  and  the  high 
cost  of  computation.  It  should  be  emphasised  that  numerical  modeling 
technology  is  such  that  it  could  predict  accurately  changes  in  movable  river¬ 
bed  profiles  only  if  sufficiently  accurate  input  data  were  given. 

Based  on  the  present  study,  several  recomnendat ions  concerning 
the  applicability  of  the  tested  models  to  the  GREAT-II  study  reach  can  be 
made.  First,  the  usage  of  CIIAR2  is  not  practical  because  the  program  is 
not  readily  available  (a  contract  with  the  Sogreah  is  necessary) .  However, 
the  CHARE  model  is  attractive  because  of  the  simplicity  of  its  model 
construction,  requiring  the  least  input  information  among  those  models 
tested,  and  its  economic  advantages.  Second,  the  CSU  models  should  be 
used  in  predicting  local  shoaling  areas  because  of  their  capability  to  pre¬ 
dict  lateral  changes  in  river-bed  elevations  which  neither  CIIAP2  nor  HEC-6 
is  able  to  provide.  Third,  I3EC-6  should  also  be  used  for  its  capability  of 
analyzing  river-bed  armoring  processes,  although  it  cannot  predict  lateral 
changes  of  bed  elevations.  Finally,  it  is  strongly  recommended  that  these 
models  be  retested  for  Pool  20  upon  compilation  of  a  new  set  of  detailed 
cross-section  data. 
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Appendix  B  Sample  HEC-6  Input  and  Output  Formats 

Input 

t i  i^wa  institute:  of  hydraulic  pcsfarch  <G323>  i i hr o  01  page: 
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Appendix  B  Cont’d 
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T 1  IOWA  INSTITUTE  OF  HYDRAULIC  RFSFAKCH  (G373)  IIHR001  M  G>'  CD. 
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468 

1450 

462 

1500 

457 

2100 

460 

22  CO 

458 

27  5  0 

GR 

4  59 

2300 

480 

2400 

M 

354.0 

50 

2350 

X 1354.54 

17 

0 

3630 

2851 

2851 

2851 

OF 

4  30.0 

0 

471.6 

26  q 

48  0.0 

515 

48  4.0 

615 

4  0  4  .  0 

765 

r.  f? 

4  P  0  •  0 

865 

477.3 

515 

4  9  0.0 

96  5 

4  0  3.0 

1  C  9  C 

400.0 

12  15 

CR 

4  77.3 

!  !  6  2 

4  75.8 

1*58 

470.6 

2  03  6 

464.9 

2035 

460.5 

31  74 

G  D 

455.3 

3515 

49  0 . 0 

3630 

h 

354.54 

1215 

3630 

X135^.ci 

13 

0 

36  0  7 

2059 

2059 

2059 

CP 

4f  3.0 

0 

472.6 

2  69 

48C.C 

5  15 

480.0 

615 

4  76.4 

7  98 

Z  r 

4  6  7.9 

1038 

464.3 

2045 

465.2 

2412 

468.4 

2773 

4  70  .4 

31  44 

r  r 

^0.0 

3318 

470.9 

3  4  7  6 

480.0 

360  7 

)* 

384. <*3 

615 

3313 

*1355.39 

15 

C 

3640 

2429 

2429 

2429 

o 

420.0 

0 

4  76.7 

126 

4  73.2 

39  0 

4  71.9 

68  2 

468.0 

?  42 

c  c 

466.0 

1  c  6  5 

4  6  4.5 

1044 

465.5 

2  39  8 

47  ^  .  0 

2677 

480.0 

7792 

GP 

4-3.0 

284  2 

4  R  %  .  0 

3292 

4  9  0.0 

334  2 

4  6  G  •  2 

35  0  0 

4  0  C  .  0 

36  31 

M 

355.39 

0 

2792 

*  1  3  5  5  .  5  9 

17 

0 

3380 

1056 

1C56 

1  OL  6 

G  P 

480 .0 

0 

470.4 

63 

463.1 

26  6 

46  6*4 

791 

465.1 

1024 

OP 

4  66.1 

1519 

4  65.8 

1700 

4  69.6 

1  99  0 

4  7  c  .  2 

22 1  P 

40C  .0 

72  97 

C' 

4  6  9.5 

24  2  5 

4  7  4.4 

26  9  3 
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269  1 

4  0  .9 . 0 

289  1 

40  0.0 

30  91 

C  R 

4  75.1 

3249 

480.0 

3380 

R 

355.59 

0 

2297 

X  1 

356.0 

7 

0 

1800 

2165 

2165 

2165 

GP 
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0 

478 

50 

460 

2C  0 

459 

350 

478 

1350 

G  P 

4  70 

1750 

43  0 

1800 

H 

356.0 

50 

1700 

X  1 

357.0 

18 

0 

7.  C  0  0 

5280 

5280 

8280 

GR 

480 

0 

4  74 

50 

472 

200 

478 

350 

473 

450 

GR 

477 

600 

4  74 

900 

477 

1  000 

466 
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20  00 

GR 

468 

2C50 

467 

2200 

465 

2250 

467 

24  00 

466 

2450 

GP 

467 

2750 

473 

250  0 

400 

300  0 

H 

357.0 

50 

2950 

X  1 

358.0 

9 

0 

1300 

5280 

5280 

5200 

GR 

479 

0 

471 

200 

473 

25  0 

465 

350 

464 

4  50 

GP 

459 

750 

461 

90  0 

457 

1250 

471 

1300 

H 

350.0 

50 

1300 

X  1 

359.0 

10 

0 

2300 

5280 

528  0 

5280 

GP 

480 

0 

471 

300 

474 

450 

472 

1200 

470 

1250 

u  R 

472 

1350 

465 

1700 

456 

1  ?5  0 

469 

2200 

480 

23  00 

H 

359.0 

100 

2250 

X  1 

360.0 

9 

0 

130  0 

5280 

5280 

^280 

GR 
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0 

460 
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45  4 
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GP 

460 

1100 

458 

1150 

473 

1200 

480 

1300 

H 

360.0 

50 

1  250 

X  1 

361.0 

1  0 

0 

1800 

5280 

5280 

5280 

CR 

481 

0 

476 

150 

470 

200 

479 

350 

471 

500 

OR 

458 

1100 

467 

1650 

466 

1700 

471 

1750 

481 

13  00 

H 

G  T 

X  1 

361.0 

400 

1800 

362.0 

10 

0 

2900 

5280 

5280 

5280 

OR 

43  0 

0 

471 

12  0  0 

4  72 

1500 

471 

1750 

475 

18  0  0 

OR 

467 

1950 

470 

2150 

470 

2650 

46  7 

2750 

480 

29  00 

F 

362.0 

1200 

2980 

X  1 

363.0 

8 

0 

2400 

5280 

“280 

5260 

GR 

480 

0 

471 

850 

470 

1  COO 

469 

1200 

468 

1500 

GP 

470 

1950 

4  /  2 

2200 

480 

2400 

H 

363.0 

t  )  0 

2400 

XI 

364.2 

14 

0 

2100 

6336 

6  33  6 

6336 

GR 

48  0 

0 

474 

950 

474 

1  050 

473 

1150 

4  74 

1250 

0  R 

473 

1350 

4  72 

1450 

466 

155  0 

466 

1600 

4  7  D 

17  00 

G  r 

469 

1850 

47  3 

1950 

475 

2  05  0 

490 

2100 

H 

E  J 
T4 

364.2 
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^OPri  INCLUDES  CLAY,  SILT. 

S  AND  , 

AN 

D  GRAVEL.  TOTAL 

SET  I  ME  NT 

LOAD 

INPUT 

T  5 

A  T 

L  •'.DI¬ 

INCLUDES 

ONLY  CLAY,  SILT 

♦ 

AND  VFS 

.  HUT,  Is 

fS  MOINES 

RIVER 

SED. 

T  6 

I  \  P 

LI  7  HAS 

CLAY*  SILT,  VFS, 

F  S  ,  ”  S 

« 

C  S  ,  VCS 

,  VFG,  A 

NO  FG. 

T  7 

T  r  F 

F  ALFT  I 

FP~MUL  A 

IS  USED. 

MODEL 

UTILIZES 

WEEKLY 

AVERAGED 

D ISCH A 

RGES 

TP 

I  1 

12 

F  OF 

W  A  T  E  R 

AND  SEDIMENT. 

CLAY 

13 

SILT 

3 

3 

I  4 

SAND 

7 

L 

G 

1  G  0  0  0 

50000 

1  5  0  C  0  0 

2CC000 

L 

CLAY 

3  5 

1  79  0 

2  6  0  8  0 

5262  0 

L 

SILT 
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5  1 5  G 

7H  180 

151694 

L 

VFS 

0.83 

42 

611 

1233 

L 

FS 

0.0 

C.O 

C.O 

0.0 

L 

MS 

0.0 

0.0 

0.0 

0.0 

L 

cs 

0.0 

0.0 

0.0 

0.0 

L 

vcs 

0.0 

0.0 

0.0 

0.0 

L 

VFG 

0.0 

C.O 

0.0 

0.0 

L 

FG 

0.0 

0.0 

0.0 

0.0 

\ 

343.2 

0.052 

G.  026 

0.996 

0.005 

0.005 

0.004 

0  .0  75 

N 

0  .475 

0.303 

0.079 

0.034 

0.016 

3  44.0 

0.052 

0.026 

0.996 

0.005 

0.005 

0,004 

0.3  75 

pj 

0  .475 

0.303 

0.075 

0.034 

0.016 

r; 

345.0 

0.052 

0.026 

0.996 

0.005 

0.005 

0.004 

C.C  75 

rj 

0  .475 

0.303 

0  .C79 

0.034 

0.016 

n 

346.0 

0.052 

C.026 

0.996 

0.005 

0.005 

0 .0  04 

0.3  75 

ft 

0.475 

0.303 

0.079 

0 .034 

0.016 

•j 

347.0 

0.052 

0.026 

0.996 

0.  C05 

0.005 

0.004 

0  .0  76 

N 

0  .475 

0.303 

0.079 

0.034 

0 .016 

348.0 

C  .052 

0.026 

0  •  v  9  9 

0.001 

C  .0  C? 

0  .0  03 

C  .1  01 

f! 

0.477 

0.26° 

0.085 

0.02° 

0.C10 

* 

348.96 

G.052 

0.026 

0  .°99 

0.001 

0.00? 

C  .  C  0  3 

P  .1  01 

r 

0  .477 

0 .288 

0.089 

0  •  0  2  r< 

0.010 

fj 

349.29 

0.052 

0.026 

0.‘?5 

0  .  o  o  G 

0.000 

c  •  r  v  i 

4  .0  48 

0.409 

0.36  1 

0.12° 

C.m 

0.01? 

p 

349.45 

0 . 05  2 

0.026 

0.978 

0.00  0 

n.oo  i 

o.  no  l 

0.0  M 

0.444 

G  .  2  7  9 

0  -  C  C3 1 

0 .068 

f.  .  r'  *  5 

N 

349.82 

0 . 0  *'  2 

0  •  r.  2  f » 

0.04  7 

G.047 

0  .  G  34 

0.1  66 

N 

0 .331 

0.19? 

0.084 

0.0*9 
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ct;c%  fed  CHfe»r.£  u*  tlcv  thal.cg  a  srriKEf.r  load  r.  to,s/c/.y 

:  Ff  £  T  F  f  t  T  EL  FEET  CF<  CLAY  MLT  £A:.P 

?64.?C0  -.05  4°S.71  465.95  Sf .**."■  7 .  9  r  79  .  2r:\  ?*■  .  r53 

262. 0C5  -.03  485.37  467.57  964*7.  ,'e7!/.  25576.  1015, 

362.0  00  -.02  4^5.17  466. r«8  °64;?7.  ►.-7«a.  25376.  1 4 1 9 

361.0  00  -.05  <*f?4.?6  457. >5  1127*W.  14-53.  4795'.  1231? 

360  . 0  C0  .06  494.6?  454.0b  1  1275  7.  14453.  47553.  11619 
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